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foreword 


This quarterly review of reactor development has been prepared at the request 
of the Division of Information Services of the U. S. Atomic Energy Commission. 
Its purpose is to assist interested organizations in the task of keeping abreast 
of new results in reactor technology for civilian application. 

The report is a concise discussion of selected phases of research and de- 
velopment for which there have been significant advances or a heightened 
interest in the past few months. It is not meant to be a comprehensive abstract 
of all material published during the quarter, nor is it meant to be a treatise on 
any part of the subject. The intention is to cover the various areas of reactor 
development from the general viewpoint of the reactor designer rather than 
from the more detailed points of view of specialists in the individual areas. 
However, papers which are thought to beof particular significance or particular 
usefulness in specialized fields will be mentioned in short notes. In the over- 
all plan of the report, it is intended that various subjects will be treated from 
time to time and will be brought up to date at that time. 

Any interpretation of results which is given represents only the opinion of 
the editors of the report, who are General Nuclear Engineering Corporation 
personnel. Readers are urged to consult the original references in order to 
obtain all the background of the work reported and to obtain the interpretation 
of the results given by the original authors. 


W. H. ZINN 
General Nuclear Engineering Corporation 








contents 








Foreword 


GENERAL RESEARCH AND DEVELOPMENT 


II 


III 


IV 


REACTCR APPLICATIONS: PROCESS STEAM GENERATION 


Uses of Process Heat 

Process Heat Generation and Consumption 
Reactors for Producing Process Steam 
References 


THERMIONIC ENERGY CONVERSION 


Basis of Operation 

Efficiency Equation 

Estimates of Attainable Efficiencies 

High-vacuum Thermionic Converter Performance 
Experiments on a Cesium-vapor-filled Thermionic Converter 
Applications 

Los Alamos Plasma Thermocouple 

References 


REACTOR PHYSICS: H,O-MODERATED CRITICAL ASSEMBLIES 
AND ANALYSES 


References 


HEAT TRANSFER AND FLUID FLOW 


Collapse of Parallel-plate Fuel Assemblies 
Nucleate Boiling 

Liquid-metal Heat Transfer 

Liquid Metals with Internal Heat Generation 
Water Flowing Parallel to Tube Bundles 
References 


REACTOR KINETICS AND DYNAMICS 


Delayed Neutrons in Reactor Kinetics 
Reduced Delayed-neutron Group Representations 
Approximations for Reactor Accident Calculations 











contents (continued) 


26 Oscillation Experiments on Large Heavy-water Reactors 
26 Transfer Functions of Distributed Parameter Systems 
27 References 


28 VI REACTOR SAFETY AND CONTAINMENT | 


28 AEC Reactor Safety Program 
30 Metal-Water Reactions 

33 SRE Shutdown 

33 References 


36 VII SHIELDING 


36 Review of Shielding Art 
36 Reference 


37 VIII MATERIALS 


37 Radiation Effects on Semiconductors 

37 Radiation Effects on Magnetic Materials 

38 Corrosion of Stainless Steel in a Nitrogen Atmosphere 
40 References 


41 | PROGRESS ON SPECIFIC REACTOR TYPES 


41 IX LIGHT-WATER REACTORS 

41 The Belgian BR-3 

41 | Reference 

42 X BOILING-WATER REACTORS: STABILITY 

42 Design Considerations 

a4 Experimental Boiling Water Reactor 

47 BORAX Reactors 

51 Review of ANL Experience 

53 References 

54] XI FLUID-FUEL REACTORS | 
34 Report of the Fluid Fuel Reactors Task Force | 
38 The Homogeneous Reactor Test 
60 References | 


61 XII NEW REACTOR CONCEPTS: THE FLUIDIZED-BED REACTOR 
STUDY 





63 References 





issued quarterly by the U. S. Atomic Energy Commission. Use of funds for printing this 
publication approved by the Director of the Bureau of the Budget on November 1, 1957. 


iv 





GENERAL RESEARCH AND DEVELOPMENT 





REACTOR APPLICATIONS: PROCESS STEAM GENERATION 





A recent Argonne report! on process steam 
generation directs attention to the question of 
nuclear process heat. 

A brief discussion of the possible future 
magnitude of nuclear process heat applications, 
based on references 2 and 3, was givenin Power 
Reactor Technology, 2(3): 4-5 (June 1959). The 
questions of more immediate importance to the 
reactor designer are the magnitude of the present 
process heat demand, the requirements of a 
process heat plant, andthe probability of meeting 
those requirements with nuclear power at the 
present stage of nuclear technology. A brief 
discussion of the first two questions, based on 
references 4 to 7, is given below; the Argonne 
report may be regarded as a consideration of 
the third. 


Uses of Process Heat 


In general, the possible applications of proc- 
ess heat fall into two broad categories defined 
by operating temperature: a low-temperature 
range of 400°F or less for use in manufacturing 
processes requiring low-quality steam, and a 
high-temperature range of 1500°F or more for 
use in certain specialized chemical and metal- 
lurgical operations. The latter applications are 
outside the range of existing reactor technology 
and will not be considered further here.* The 
practicability of using nuclear reactors for the 
generation of low-temperature process steam 
depends on economic factors since the technical 
feasibility of reactor operation under somewhat 
more severe conditions is well established. 

Large quantities of low-temperature process 
steam are presently consumed inthe processing 





*Some specific high-temperature applications are 
discussed in references 4 to 6. 


of products such as textile fibers, chlorine and 
alkalis, fertilizers, petroleum, antibiotics, pulp 
and paper, aluminum, and corn products. In the 
majority of applications, the amount of heat 
required for an individual plant is under 80 Mw(t), 
and the user pays less than 70 cents per 10° Btu 
for energy in the form of steam. Although the 
low unit capacities of the plants do not favor the 
use of nuclear power, there are certain charac- 
teristics of these applications which are favor- 
able. Almost all the plants operate at a load 
factor of 90 per cent or greater, which is, of 
course, desirable for a nuclear plant because 
of the comparatively high contribution of capital 
investment to the total steam production cost. 
A second consideration is that the low tempera- 
ture and pressure requirements may allow 
cheapening of the reactor. They may, for 
example, make possible the use of inexpensive 
materials, such as aluminum or its alloys inthe 
fabrication of the core, or of carbon steel in- 
stead of stainless steel in the primary system. 
Finally, many of the industries using process 
steam are located in remote areas near their 
sources of raw materials. The use of a nuclear 
reactor would provide greater flexibility in the 
location of plants since transportation of fuel 
would not be a great problem. The remote 
location might also conceivably ease the re- 
quirements on containment. 


Process Heat Generation 
and Consumption 


Certain other aspects of the process steam 
application appear unfavorable to the introduc- 
tion of a new energy source. Since in most of 
the applications proposed the cost of the steam 
is a small fraction of the cost of producing the 
end product, proven reliability of the steam 


1 








generator is usually mandatory. If standby equip- 
ment must be provided to ensure against loss of 
production in case of reactor failure, the eco- 
nomic picture probably becomes quite dis- 
couraging. Furthermore, the potential cus- 
tomers for process steam have no great incentive 
to pursue or support any developmental pro- 
grams for steam generators, since their in- 
terests lie primarily in their finished products. 

Reference 7 gives more detailed data on the 
quantity of process steam used in past years, 
its cost, and unit sizes of installations. In 1954, 
approximately 8.7 x 10'° Btu of energy was used 
as process heat in manufacturing industries in 
the United States. Of this, about one-third was 
used in the application of direct firing of fuels 
to the process, that is, in general, at high tem- 
peratures. The remaining two-thirds was used 
as steam. The process energy used as steam 
corresponds to approximately 180,000 Mw-years 
per year, an amount roughly equal to the energy 
consumption for the production of electricity in 
that year. Most of the installations fell in the 
size range below 250,000 lb of steam per hour 
(~80 Mw(t)), and only about 5 per cent of the 
capacity was in units having individual capacities 
greater than 500,000 lb per hour (~160 Mwi(t)). 
The distribution of boiler capacity for process 
steam units added from 1945 to 1956 is given, 
over the range of steam pressures and unit 
capacities, in Table I-1. 


Table I-l1 PERCENTAGES OF BOILER CAPACITY ADDED 
IN VARIOUS CAPACITY AND PRESSURE RANGES’ 
(United States Manufacturing Industries, 1945 to 1956) 





Percentages of boiler capacity added for the 
steam capacity ranges (in 10° lb/hr) indicated 








Pressure 

range, psi Under 100 100-250 250-500 500+ Total 

Under 250 29.1 5.5 34.6 

250-600 8.8 11.8 3.6 0.5 24.7 

600-1000 2.3 19.2 10.1 0.2 31.8 

1000 + 0.1 0.4 4.2 4.2 8.9 
Total 40.3 36.9 ato 4.9 100.0 





It is estimated’ that, for the smaller coal- 
fired process steam boilers operating at 90 per 
cent load factors, the cost of steam generation 
is about 70 cents per 10001b, whereas for a large 
boiler in the 600,000 lb/hr capacity range, the 
generating cost under the same conditions 
would be about 47 cents per 1000 lb. These costs, 
in terms of energy, correspond to approximately 
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65 cents and 44 cents per million Btu, or 2.2 and 
1.5 mills/kw(t)-hr, respectively. 


Reactors for Producing 
Process Steam 


The Argonne report! covers a feasibility study 
of nuclear reactors for process steam which 
was undertaken at the request of the Atomic 
Energy Commission. The area defined for study 
was characterized by the production of saturated 
steam at 380°F and ata rate of 125,000 lb/hr, 
with shielding and containment of the plant suit- 
able for location on an industrial site. Four 
reactor types were studied: a pressurized-water 
reactor (PWR) with stainless-steel primary 
system, a similar PWR with carbon-steel pri- 
mary system, a boiling-water reactor (EWR), 
and an organic-moderated reactor (OMR). The 
estimated steam production costs (in cents per 
1000 lb of steam) for the designs studied were as 
follows: 


‘PWR 64-150 
BWR 83-172 
OMR 75-154 


The PWR listed is the one employing the 
stainless-steel primary system. The ranges in 
production costs result primarily from differ- 
ences in assumptions as to rates of capital 
charge and plant factors. Table I-2 gives a 
breakdown of production costs for a capita] 
charge rate of 15 per cent per year and a plant 
factor of 90 per cent. It is stated in the report 
that none of the reactor types studied has 
dominant characteristics which show it to be a 
clear-cut choice over the others for the appli- 
cation. It is implied that the considerably higher 
fuel cost for the boiling-water reactor is a 
consequence of the specific reactor design 
chosen rather than a characteristic of the 
reactor type. 

Other conclusions reached in the study are that 
carbon or low alloy steels should be acceptable 
materials of construction for water-cooled reac- 
tors having operating temperatures of 500°F or 
lower, provided that the pH is maintained at 10 
to 11, that hydrogen is added to the water, and 
that dissolved and suspended solids are main- 
tained at low levels. It is also concluded that 
aluminum-clad fuels with surface temperatures 
as high as 500°F maybe acceptable in stainless- 
steel systems where the fH of the water can be 
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Table I-2_ STEAM GENERATION COSTS! 

















PWR, 
stainless PWR, 
steel carbon steel BWR OMR 
Capital cost of plant, including land and escalation, 
dollars 2,763,000 2,486,000 2,226,000 2,260,000 
Capital charge at 15 per cent per year and 90 per 
cent plant factor, cents/10° lb steam 42 38 34 34 
Nonfuel operating charge, cents/10° lb steam 15 15 15 21 
Fuel operating charge, cents/10° lb steam 20 20 28 24 
Fuel inventory charge, at 4 per cent per year on 
uranium, 6 per cent per year on fabrication cost 
of fuel inventory not yet loaded, and 15 per cent 
per year on fabrication cost of core in reactor 12 12 32 17 
Total, cents/10° lb steam 89 85 109 96 
Table I-3 SUMMARY OF PROCESS-STEAM REACTOR DESIGN DATA! 
PWR BWR OMR 
Power, Mwi(t) 40.6 40.6 40.6 
Fuel UO, in Zr UO, in Zr U plates in Al 
Coolant H,O H,O Diphenyl 
Core diameter x length, in. 42 x 42 56.3 x 60 48 x 48 
Core volume, liters 953 2442 1475 
Power density, kw/liter of core 42.6 16.6 27.5 
Fuel enrichment, % 1.8 1.5 is 
Initial fuel loading, kg of U 3380 10,300 7650 
Specific power, kw/kg of U 12.0 3.9 5.3 
Burn-up, Mwd/metric ton 17,200 12,300 8900 
Nonwater to water ratio in fuel region 1.35:1 1.231 
Nonorganic to organic ratio in fuel region a=% 
Initial conversion ratio >0.8 >0.8 >0.8 
Void coefficient of reactivity, 
ae. pe... ~ 0.2284 0.1884 
k water volume 
Average steam volume fraction in core 0.18 
Average heat flux, Btu/(sq ft)(hr) 9.08 x 104 4.36 x 104 1.65 x 104 
Maximum heat flux, Btu/(sq ft) (hr) 4.74 x 105 1.30 x 10° 7.6 x 104 
Primary flow rate, lb/hr 4.22 x 108 1.39 x 10° 7.9 x 10° 
Inlet temp., °F 400 239* 421 
Outlet temp., °F 430 427 450 
Reactor pressure, psi 900 333 75 
Process steam conditions: 
Flow rate, lb/hr 1.25 x 10° 1.25 x 10° 1.25 x 10° 
Temp., °F 380 380 380 
Pressure, psia 196 196 196 
Condensate return temp., °F 150 150 150 
Make-up water, lb/hr at 60°F 35,000 35,000 35,000 





*Inlet temperature into reactor vessel. 


5 to 6, but further experience with aluminum 
corrosion is needed before a definite recom- 
mendation for its use under these conditions can 


be made. 

The main characteristics of the designs ar- 
rived at for the three reactor types are given in 
Table I-3. The pressurized-water reactor uses 
a two-pass flow system through the core, and 


its heat-removal characteristics are defined by 
the condition that no fuel-element surface tem- 
perature shall be higher than the water satura- 
tion temperature. The boiling-water reactor isa 
natural-circulation system. The process steam 
is generated in a reboiler whose primary side 
is fed only by primary steam from the reactor. 
The reactor employs absorbing control rods with 








4 GENERAL RESEARCH AND DEVELOPMENT 


fuel bundles as followers, but these do not appear 
to have any great effect on the economic evalua- 
tion. The organic reactor is moderated and 
cooled by diphenyl; it employs a plate type fuel 
element of uranium —3.5 wt.% molybdenum—0.5 
wt.% silicon alloy, nickel-plated, and clad with 
finned 2S aluminum. Both of the water reactors 
are provided with steel secondary containment 
shells; the organic reactor is housed in an 
insulated sheet-steel building with reinforced- 
concrete basement (in which the reactor and its 
primary system are located) and special sealing 
to minimize leakage through the walls of the 
steel building. 

If a steam production cost of about 70 cents 
per 1000 lb can be considered a reasonable 
economic target for a reactor of this size, then 
the study would appear to indicate that the 
reactors considered have some chance of ap- 
proaching this goal. Since it is probable that 
there is a greater variation in process-steam 
generating costs than in electric power generat- 
ing costs in this country, it may well be that 
there are individual instances in which the use 
of nuclear process steam would appear eco- 
nomically attractive. This, of course, could be 
determined only by making a more detailed 
reactor design for the particularly favorable 
case selected. 

In view of the simplicity of the natural- 
circulation boiling reactor and its lowoperating 
pressure for the process-steam application, and 
in view of the encouraging experience that has 
been had with the construction costs and operat- 
ing behavior of EBWR and other experimental 
boiling reactors, it comes as a surprise that 
the study should show a substantially higher 
steam production cost for the boiling reactor 
than for either the pressurized-water or the 
organic-moderated reactor. Although this dif- 
ference is de-emphasized in the report as 
being attributable to the specific design study, 
it must nevertheless be presumed to have some 
meaning if the study has significance as a 
comparative study. A consideration of the cost 
breakdown of Table I-2 shows that the boiling 
reactor does have the lowest estimated capital 
cost, that most of the excess cost of steam 
generated by the boiling reactor is due toa very 
high fuel inventory charge, and that the re- 
mainder of the excess is due toahigh fuel oper- 
ating charge. The latter difference between the 
boiling and the pressurized reactors appears to 
stem from the fact that the two reactors were 


designed with equal quantities of initial excess 
reactivity, with the fuel lifetime defined as the 
reactivity lifetime in a uniform flux (obtained 
by radial and axial fuel shuffling), single-batch, 
reloading cycle. This appears to be a rather 
arbitrary way of defining the fuel lifetime, un- 
less it could be shown that a real limitation is 
involved which could not be circumvented either 
by the provision of more control-rod capacity or 
by the use of a partial reloading scheme. 

The very high inventory charge against the 
boiling reactor reflects the low specific power 
of 3.9 Mw(t) per metric ton of uranium (Table 
I-3), which in turn is consistent with the low 
power density of 16.6 kw/liter of core. These 
low values evidently result from the high value 
of fuel-to-water ratio used in the reactor. The 
ratio of nonwater to water volume in the core is 
1.2. Thus the power density per liter of coolant 
is 16.6 x 2.2 =37 kw/liter, a value which is not 
high for a boiling reactor but nevertheless is not 
so strikingly low as the power density per liter 
of core. The second effect of the high fuel-to- 
water ratio is a relatively high negative steam 
void coefficient of reactivity (-0.188 per cent 
keg¢/per cent void). As a result of this rather 
high coefficient, the reactivity held by steam 
apparently amounts to 3 per cent Reff or more 
even at the low design power density. The high 
fuel-to-water ratio was evidently used in order 
to attain a high conversion ratio. Indeed, the 
conversion ratios quoted in the report for the 
three reactor types— which were computed as 
part of a parametric study rather than for the 
specific reactor designs described—are sur- 
prisingly high for reactors of such small size. 

The conflicting demands of conversion ratio 
and steam void coefficient of reactivity have 
been pointed out before in this Review, Power 
Reactor Technology, 2(1): 16 (December 1958) 
and 2(4): 13-16 (September 1959), and are dis- 
cussed further in the section on boiling reactor 
stability in the present Review, p. 41. Although 
a compromise must usually be made between 
conversion ratio and power density, and the 
optimum may very well vary with the specific 
reactor design and application, it would seem 
improbable that the optimum would swing so far 
toward high conversion ratio as to result in an 
inventory charge which accounts for some 30 
per cent of the energy generation costs. If the 
boiling reactor in question were designed for a 
“wetter” core and were given Sufficient excess 
reactivity to achieve a fuel life equal to that in 
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the pressurized-water reactor, the inventory 
charge would decrease, the net fuel burn-up 
cost would increase to some extent, and there 
should be some decrease in capital cost. The 
net effect would probably be a reduction of steam 
production cost. 
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THERMIONIC ENERGY CONVERSION 





Basis of Operation 


The thermionic energy converter is a device 
which converts heat energy directly into elec- 
trical energy. It is shown in its simplest form 
in Fig. 1. Basically, the thermionic energy 
converter consists of a hot cathode which emits 
electrons thermionically and a cold anode which 
collects the emitted electrons. The energy in- 
put to the device shown in Fig. 1 is heat which 
is directed into the cathode. From this heat, 
useful electrical power can be obtained, under 
the proper circumstances, by connecting a 
resistive load between the cathode and the 
anode. The operation of the converter can be 
viewed as a thermodynamic cycle wherein the 
working fluid is the electron gas which is made 
to flow from the hot plate to the cold plate with 
a resultant electrical work output. 


HOT CATHODE COOL ANODE 






































+ LOAD _ 














Figure 1—The simplest form of a thermionic con- 
verter. 


The thermionic converter utilizes the so- 
called “Edison effect,’’ which is probably the 
oldest observed electronic phenomenon. It has 
been only in the past two years, however, that 
the device has received serious attention as a 
possible efficient and practical heat engine. 


The basis for operation of the thermionic 
converter may be better understood by refer- 
ring to the potential-energy diagram of Fig. 2. 
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Figure 2—A potential-energy diagram ofa thermionic 
converter. 


The vertical distance on the diagram represents 
the potential energy of the electron. Heat in- 
creases the temperature of the electrons in the 
cathode, and, therefore, some of the electrons 
acquire sufficient energy to escape from the 
plate. The minimum energy required for es- 
cape is ¢-, which is the energy difference be- 
tween the Fermi energy level of the electron 
within the plate and the energy of the electron 
at rest in a vacuum outside the plate. An ad- 
ditional amount of energy must be supplied to 
the emitted electron to hurdle the space-charge 
potential barrier, ¢,. Hence, at position a in 
Fig. 2, the electron has potential energy (rela- 
tive to the cathode Fermi energy level) equal 
to dc + dp. In addition, the electron has ki- 
netic energy 2kT;, which is only a few per 
cent of the potential energy. In effect, then, 
heat has lifted the electrons over the cathode 
work-function barrier and the space-charge 
barrier to the high potential required for the 
electrons to cross over to the anode. If the 
electron falls back into the cathode, this energy 
is given back to the cathode. However, from 
position ain Fig. 2, the electrons can slide 
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down the potential energy curve toward the 
anode, thereby converting potential @p into ki- 
netic energy of the electrons. This kinetic 
energy goes into the heating of the anode when 
the electrons strike the anode. At the anode 
surface, the electrons fall into the metal and, 
as a result, lose potential @, as the electrons 
join those at the Fermi energy level of the 
anode. Note that ¢4 is the work function of the 
anode. At the Fermi energy level of the anode, 
the electrons have potential energy relative to 
the Fermi energy level of the cathode; it is 
this potential difference V), which is available 
as electrical voltage to do work in an external 
circuit. 


The desirable arrangement is a high cathode 
work function ¢, and a low anode work function 
o,4- However, if the two plates were at the 
same temperature, the anode would then emit 
more electrons than the cathode. Hence the 
anode must be kept colder. With plate spacings 
of the order of 0.01 in. in a vacuum, the space- 
charge barrier limits the electron flow to a 
very small value. In this case, ¢, is so large 
that, at practical cathode temperatures, only a 
small fraction of the electrons emitted have 
sufficient energy to penetrate the space-charge 
barrier. As the space-charge barrier is re- 
duced, the electron flow is increased and is 
limited finally by the emission capabilities of 
the cathode at its operating temperature. One 
way of reducing the space-charge barrier is to 
place the anode and the cathode very close to- 
gether.’ Another method is to introduce heavy 
positive ions to neutralize the electron space 
charge.’ Introduction of a positively charged 
grid between the cathode and the anode would 
also be helpful if the grid did not draw exces- 
sive current. Electric and magnetic fields can 
be employed to conduct the electrons from the 
cathode to the anode as described by Hat- 
sopoulous.° The space-charge barrier prob- 
lem is one major obstacle which must be 
overcome in order to make the thermionic 
device look promising for further development. 


Efficiency Equation 


This Review is based mainly on the infor- 
mation given in references 1, 4, 6, and 7. The 
Los Alamos plasma thermocouple,*” which is 
also a thermionic device, is treated in a sepa- 
rate section. 


The general analysis of reference 7 estab- 
lishes approximate upper limits—which may 
be achieved in practice —on the efficiency of 
thermionic converters. The analysis also pro- 
vides a guide as to which experimental pa- 
rameters can be profitably varied to improve 
efficiency. No details of the space-charge po- 
tential barrier are considered in the analysis; 
instead, the space-charge considerations are 
circumvented by analysis of the problem in 
terms of the potentials V. and V, indicated in 
Fig. 2. 

The efficiency of a thermionic converter is 
simply the ratio of the electrical power output 
to the power input, as given by the following 
equation:””' 


(V_ = Vale —dq) rr Vy (Je -Jq) 


sear RT RT, 
Je(Ve +! - 2° Jal Ve + =| +K+R+I1 








The numerator is the product of the output 
voltage (V. — V, ) and the net current (J, —Jg), 
minus the joule losses in electrodes and leads, 
Yw(Jc —~Jqa)*. Here, J. is the electron-emission 
current from the cathode to the anode, J, is 
the back-emission current from the anode, 7, 
is the electrical resistance of electrodes and 
leads, and (V. — Vq) is as indicated in Fig. 2. 
The first term in the denominator represents 
the energy required to lift the electrons from 
within the cathode material and over the space- 
charge barrier. The second term is the energy 
returned to the cathode by anode back-emission. 
Here, e is the electronic charge, T,. and T, are 
the cathode and anode temperatures, respec- 
tively, and k is Boltzmann’s constant. Theterm 
K represents the heat-conduction losses from 
the cathode, the term A represents the radiation 
losses from the cathode, and the term / repre- 
sents the energy input required to produce 
positive ions for space-charge neutralization. 


Estimates of Attainable 
Efficiencies 


Houston’ first treats the idealized case of no 
radiation and no conduction losses and where 
there is zero energy input for positive ion pro- 
duction. His results show that there exists an 
optimum value of V, for specified V-, 7; , and 
Tq values. At this optimum, the efficiency is 
roughly 90 per cent of the Carnot efficiency. 
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Note that since the device is a heat engine, its 
efficiency can never exceed the Carnot effi- 
ciency as given by T, — T,/T,. At high anode 
temperatures (7, ~1000°K), the optimum V, 
is in the neighborhood of 1.5 volts; this is 
within the realm achievable with presently 
known materials. However, for anode tem- 
peratures in the vicinity of room temperature, 
the optimum V, value is near 0.5 volt, which 
is unattainably low. 

Houston’ also considers more practical cases 
in which the extraneous energy losses are 
included in analysis of over-all efficiency. An 
anode temperature 7, of 500°K was assumed. 
However, this assumption is unimportant in 
that the results obtained are valid for fairly 
wide variations in anode temperature from the 
assumed value (7, should be low enough so that 
back-emission is negligible). For a givenvalue 
of V, and J, (where J. is inamperes per square 
centimeter of cathode surface area), optimum 
values of T, and V; are shown to exist. If a 
Vq value of 2 volts can be achieved, amaximum 
7 of 16 per cent is indicated as possible at J, = 
10 amp/cm’. If a V, value of 1 volt were pos- 
sible, the maximum 7 would increase to 34 per 
cent atJ, = 10 amp/cm’. For V, = 2 volts, 
values of T. = 2200°K and V, =3.4 volts are 
indicated as optimum. A conclusion is that a 
low value of V, and, hence, of the anode work 
function ¢4 and space-charge potential barrier 
¢», is vital in obtaining high efficiency from a 
thermionic converter. 

It is evident that a crucial question arises in 
the attainment of the lower limit of V;. Fora 
cesium-on-oxidized-silver surface, Houston’s' 
interpretation of Koller’s? measurements indi- 
cates $4 =1 volt at 400°K. Since V, = 4+ 
$», there seems to be little hope of achieving 
aV, value less than 1 volt with presently known 
materials. Hence the value of 7 = 34 per cent 
which is calculated as maximum for an as- 
sumed V, =1 volt may represent an upper 
bound of the achievable efficiency. 


High-vacuum Thermionic 
Converter Performance 


The results of calculations of the perform- 
ance of a high-vacuum thermionic energy con- 
verter are presented in reference 1. The cal- 
culations include the effect of the space charge 
as determined from Langmuir’s work.’ The 


pertinent features of the analysis and the re- 
sults obtained are outlined below. 

To determine the performance of the con- 
verter, it is necessary to know how the output 
voltage V, varies with current / flowing through 
the device. The anode is assumed to be cold 
enough so that no electrons are emitted from 
it. Then, with reference to Fig. 2, 


Vo = (dc — ba) + (be — bp) 


Both ¢, and ¢, are functions, respectively, of 
the cathode and anode materials used and are 
independent of current flow. However, the 
space-charge potentials, ¢, and $s, vary with 
current. 

If the cathode temperature 7, , the saturated- 
emission density J], available from the cathode, 
and the separation d of the cathode and anode 
are known, then both ¢, and ¢, can be calcu- 
lated'® as functions of J. The generalized re- 
sults obtained were used to predict the per- 
formance that can be expected of energy 
converters operating at7. =1000 to 1200°K 
with saturated-emission densities character- 
istic of oxide-coated cathodes run in clean 
environments.'! Figure 3, taken from refer- 
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Figure 3 — The output power density versus the current 
density for four different cathode-anode spacings.’ 
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ence 1, is a plot of output power density (watts 
per square centimeter of cathode surface area) 
versus current density. First, there is a par- 
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Figure 4—The output power density into a matched 
load versus cathode-anode spacing for various work- 
function differences.’ 


ticular value of external load resistance for 
which maximum power is transferred to the 
load. Second, the power density delivered to 
the load increases rapidly as the cathode- 
anode separation distanced decreases. In Fig. 
4, output power density at maximum power 
transfer, i.e., using a matched load, is plotted 
versus cathode-anode spacing d for various val- 
ues of work-function difference ranging between 
0 and 1 volt. The output power density in- 
creases rapidly with increase in work-function 
difference (¢, — $4). 





Webster’ points out that some of his predic- 
tions have been checked experimentally in 
energy converters built by Beggs,'! and satis- 
factory agreement was obtained. Beggs’ work 
on construction techniques of thermionic con- 
verters makes the small cathode-anode spacings 
(less than 0.001 in.) feasible. An estimate of 
the efficiency of the high-vacuum thermionic 
converter is given by Webster.' With 7, = 
1100°K, J; = 5.0 amp/cm?, d = 0.00044 in., and 
o- — 64 =0.7 ev, Webster quotes an efficiency 
of 10 to 20 per cent as attainable with average 
emissivities and average conductive heat losses. 


Experiments on a Cesium-vapor- 
filled Thermionic Converter 


Reference 4 gives the results of experiments 
on an experimental thermionic converter which 
utilizes partially ionized cesium vapor for 
space-charge neutralization. As indicated in 
Fig. 2, the space-charge potential (¢;) con- 
stitutes a barrier to the flow of electrons from 
the cathode to the anode. If the electrons are 
moving through a vacuum, the total space 
charge is simply the sum of the negative elec- 
trical charges on all the electrons in the space 
between the cathode and the anode. It will be 
distributed in this space in a manner deter- 
mined by the variation in electron velocity 
across the space, which in turn depends upon 
the potential distribution resulting from the 
space-charge distribution. If the electron cur- 
rent across the space is increased simply by 
increasing the emission from the cathode, the 
number of electrons in the space will increase 
and the space-charge and the resulting nega- 
tive potential resulting from it will also in- 
crease. As the space-charge potential becomes 
higher and higher, only the more energetic 
electrons emitted by the cathode can cross the 
space. Thus the space-charge potential im- 
poses a limit on the current density across the 
space which may be lower than the limit on the 
density of electron emission from the cathode. 

If positive ions can be introduced into the 
space, their positive charges will tend to neu- 
tralize the negative electron space charge, and 
the current-limiting effect can be greatly re- 
duced. Since the positive ions, with their much 
greater masses, move across the space much 
slower than the electrons, a positive ion charge 
equal to the electron space charge can be 
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Table II-1 RESULTS OF EXPERIMENTS ON A CESIUM-VAPOR-FILLED THERMIONIC CONVERTER‘ 
Input power Output power 
Correction Load 
Cs Cs Ion for half- Total resist- Poten- Power 
temp., pressure, Cathode, source,* time output, input, ance, tial, Current, output, Efficiency, 
Condition a. © mm Hg watts watts watts watts ohms volts amps watts % 
1 261 0.6 5.35 7.35 1.87 14.57 1.0 1A 1A 1.2 8.3 
2 261 0.6 8.80 i 0.75 9.55 2.5 LA 0.44 0.48 5.1 
3 291 1.5 5.35 7.35 1.22 13.92 1.0 0.72 0.72 0.52 3.7 
4 291 1.5 8.80 tT 2.09 10.89 0.6 0.78 1:3 1.0 9.2 





* Estimated input power for the central 2 cm of water. 
t+ Tube operated as a diode with only the cathode hot, 


maintained with an ion production rate which 
is very much smaller than the rate of electron 
emission from the cathode. The attractive 
method of producing positive ions is by the 
introduction of neutral atoms whose ionization 
potential is low enough that they may be ionized 
by thermal energy at the temperature existing 
in the thermionic device. Among the chemical 
elements, the alkali metals are the most easily 
ionized; of these metals, cesium has the lowest 
ionization potential (3.87 volts). Consequently, 
cesium vapor is an attractive substance for 
promoting space-charge neutralization. The 
presence of the cesium also affects the work 
functions of the two electrode surfaces so that, 
by adjustment of the cesium vapor pressure, 
the work functions of the cathode and the anode 
can be varied. On the cold anode the cesium 
presents a low-work-function collector. A de- 
scription of one of the several thermionic 
converter tubes built is presented in reference 
4. It is a three-electrode tube consisting of a 
cylindrical anode and two hot wires running 
through the cylinder. One wire (of molybdenum) 
is the cathode when the tube is operated asa 
diode. For operation as a triode, the other 
wire (of tungsten) is run hotter than the mo- 
lybdenum wire and serves as an ion source. 
The anode surface is a cesiated silver oxide 
surface, described by Koller,’ which is the 
lowest work-function photoemissive surface. 
The tube is enclosed in an electric oven main- 
tained at an adjustable elevated temperature to 
produce the desired cesium pressure. The 
anode is cooled by radiation and by cesium 
convection to the glass wall of the tube. 

Table II-i1 shows the results of some tests 
on this tube. It is evident that the power con- 
sumed by the ion source accounts for a sub- 
stantial fraction of the input power in those 


cases when the ion source is used (conditions 
1 and 3). If, for example, in condition 1 the 
power to the ion source is not included in com- 
puting the efficiency, the efficiency increases 
to 16.8 per cent. This suggests what can be 
gained by a more efficient ion generator or by 
making more efficient use of the ions that are 
formed. For condition 4 in the table, the open- 
circuit voltage is 2.75 volts. Hence the internal 
voltage drop is about 2.75-0.78 =1.97 volts 
and suggests that the power dissipated in the 
plasma is 1.97 x 1.3 = 2.6 watts, which is 2.6 
times the output power. Any reduction in the 
plasma potential will appear as an increase in 
the output potential and in the over-all effi- 
ciency. If, for example, the plasma drop could 
be reduced to 0.2 volt, the efficiency would in- 
crease from the measured value of 9.2 to 3C 
per cent. Unfortunately, this appears to be 
difficult to achieve. At present, it is impos- 
sible to give an accurate physical picture of 
low-voltage plasma. The true source of the 
positive ions, as well as their lifetimes and 
their concentrations in a temperature and volt- 
age potential field, is unknown. 


Applications 


Possible applications of the thermionic energy 
converter are discussed in reference 4. The 
thermionic converter is a high-temperature 
heat engine whose optimum cathode tempera- 
ture appears to be between 1300 and 2500°K. 
The anode temperature must be considerably 
lower, but there appears to be no advantage to 
be gained by anode temperatures below 650°K. 
Since present coolants in nuclear reactors 
operate in the 600°K range, the thermionic 
converter could conceivably operate as a top- 
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ping device for nuclear reactors. One way 
this might be accomplished is through the use 
of a reactor fuel element having a double- 
walled jacket with cesium vapor in a gap be- 
tween the two walls. The inner wall would be 
the cathode heated by fission, and the outer 
wall would be the anode cooled by the conven- 
tional coolant. Direct production of electrical 
energy would accompany the flow of the fission 
heat through the cesium-vapor space to the 
coolant. The heat absorbed by the coolant 
would be converted to electricity by means of 
the conventional heat cycle withthe conventional 
coolant and turbine-generator system. 

As an example of the potentialities of such a 
topping device, one might consider athermionic 
converter efficiency of 20 per cent, a value 
about halfway between that presently attained 
in working models and that indicated as theo- 
retically possible. All of the reactor fission 
heat flows through the device, and, hence, 20 
per cent of this heat appears as direct elec- 
trical energy. Eighty per cent of the heat is 
rejected to the coolant. If this heat is utilized 
with 30 per cent efficiency in a conventional 
system, the resulting electrical energy amounts 
to 0.8 x 30 = 24 per cent of the original fission 
heat. Hence the total efficiency of the system 
would be 20 + 24=44 per cent. The practical 
problems involved in using thermionic conver- 
sion in the reactor are, of course, numerous. 

Other applications might result from the 
need for electrical power sources with no 
moving parts either for very long life or for 
quiet operation. Reference 6 presents a pre- 
liminary evaluation of the thermionic power 
plant with a nuclear reactor as heat source for 
outer space applications. The heat is trans- 
ferred by conduction from the reactor core to 
the reactor surface and thence to the cathode 
surface by radiation. The reactor must there- 
fore be capable of high-temperature operation 
(4000 and 3000°F as temperatures of reactor 
centerline and surface, respectively). The an- 
ode is cooled by radiation to the space environ- 
ment. Weight considerations favor relatively 
small radiator areas, and, consequently, high 
anode temperatures are required (2200°F in 
the reference design). The electrical power 
output of the power-plant package is limited to 
about 30 kw(e) by the conduction heat-transfer 
process from reactor core volume to reactor 
surface. The efficiency of the thermionic con- 
verter power plant was estimated as 9.7 per 





cent. The specific power in electrical watt- 
output per pound of total power-plant weight 
was estimated as 18.8 (not including weight of 
any shadow shielding that may be required for 
an instrumented space probe). 

The thermionic converter has interesting 
potential applications as a simple small elec- 
tric generator converting heat from a high- 
temperature flame. In this case its operation 
is independent of an electrical distribution 
system. For example, such a device could use 
the exhaust flame in an airplane or missile for 
auxiliary power generation. In some such ap- 
plications the energy storage problem becomes 
important. 


Los Alamos Plasma 
Thermocouple 


A program has been under way at the Los 
Alamos Scientific Laboratory to develop a 
thermionic energy device for converting heat 
from fissioning material directly into elec- 
tricity.'*:> The converter is referred to as a 
“plasma thermocouple.’’ From a general point 
of view, the thermionic cell can be considered 
as a special case of thermoelectric converter 
with the cesium plasma forming one branch of 
the thermocouple. The hot junction is the inter- 
face between the cathode and the cesium gas, 
and the cold junction is the interface between 
the anode and the cesium gas. This alternate 
representation of the energy converter is given 
in Fig. 5 (from reference 14). 

Two experimental thermionic cells have been 
constructed and operated. In the first, indirect 
heating of the tantalum emitter plate is ac- 
complished by electron bombardment with an 
electron-gun beam. In the second, the emitter 
is a ZrC:UC rod within which fission heat is 
produced by insertion in the core of the Omega 
West Reactor (OWR) at Los Alamos. A cross 
section of the fission-heated thermionic cell 
is given in Fig. 6. Reference 12 presents a 
short review of the initial work and results 
obtained on the Los Alamos plasma thermo- 
couple. Since the reference is readily avail- 
able, the information given therein is not 
reviewed here. 

Reference 15 presents data obtained on the 
vacuum emission characteristics of UC, ZrC, 
and (ZrC)o.3(UC)o.2., Uranium carbide is found 
to be an exceptionally good emitter over the 
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Figure 5—Representation of a thermionic device as 
a thermocouple. 7; = hot junction at cathode tempera- 
ture, T, = cold junction at anode temperature. 
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Figure 6— Fission-powered thermionic device oper- 
ated at Los Alamos. 


observed temperature range from 1400 to 
1800°K. Pure UC is, however, unsatisfactory 
as a structural element because of its poor 
thermal-shock resistance. By contrast, ZrC 
is much more resistant to thermal shock and 
thermal cycling. The same resistance to shock 
and cycling is observed for ZrC:UC mixtures. 
The experiments show that a solid-solution 
form of uranium carbide and zirconium carbide 
gives substantially the same thermionic emis- 
sion as uranium carbide. Compatibility of UC 
and ZrC with cesium vapor was said to be 
demonstrated, although no details are included 
in reference 15. At 1800°K, the electron emis- 
sion for ZrC:UC is 1600 times higher than for 
tantalum (in a vacuum). 

Short-circuit current, obtained by placing a 
low-impedance ammeter directly across the 
terminals of a thermionic cell, is a useful guide 
to cell performance. Roughly speaking, the 
current delivered to a matched external load 
is about one-half the short-circuit current. 
Reference 14 presents short-circuit current 
data for cells containing three different emitter 
substances: Ta, ZrC, and ZrC:UC for a range 
of emitter temperatures from 1600 to 2600°K 
and of cesium pressure from 10°’ mm Hg to 
2 mm Hg. At the lowest pressure to 10°° mm 
Hg, the short-circuit current with the tantalum 
emitter is almost identical with the saturated 
emission from the tantalum surface in vacuum. 
As the cesium pressure is increased, the 
cesium atmosphere greatly enhances the emis- 
sion of tantalum. At 4x 10°’ mm Hg, short- 
circuit current exceeds saturated emission by 
a factor of 100 at 1800°K and by a factor of 
3 at 2600°K. 

No enhancement of current emission by the 
cesium atmosphere (at 3 x 107! mm Hg) occurs 
with the ZrC:UC emitter. Up to 2000°K, the 
short-circuit data agree with the saturated- 
emission data, as shown in Fig. 7. The solid 
curve represents the measured Ssaturated- 
emitter curve in vacuum; the dashed curve is 
an extrapolation of the vacuum emission data; 
the data points represent the measured short- 
circuit densities at 3x10! mm Hg cesium 
vapor pressure. At 2000°K, the ZrC:UC cell 
produces 40 times the current density of the 
tantalum cell operated under the same con- 
ditions. On this basis it is concluded that the 
plasma drop is negligible for the current den- 
sities obtained in the tantalum cell at 2000°K 
and, hence, that the tantalum-cell current at 
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Figure 7—-Comparison of saturated current emission 
as measured in vacuum (solid curve and extrapolated 
dashed curve) with short-circuit current for cell with 
ZrC: UC emitter operated with cesium atmosphere. 


2000°K is wholly emission limited. Also, from 
the agreement between saturated-emission and 
short-circuit current at values as high as 
3 amp/cm’, it is tempting to say that even in 
the ZrC:UC cell the current is still primarily 
emission limited. 

At 2100°K and above the short-circuit cur- 
rent is substantially depressed below that ex- 
pected from the saturated-emission curve. 
Reference 14 considers this as a first direct 
demonstration of a plasma-limited current 
rather than an emission-limited current. 

The highest short-circuit current observed 
at 2620°K is 124 amp for the 2-cm’ emitter 
area. The power density is about 30 watts/ 
cm’. The efficiency of the device is estimated 
as 15 per cent. 

Reference 16 discusses the thermoelectric 
properties of a gas-filled diode containing a 





high density of positive ions. The thermoelec- 
tric power of the metal with respect to the 
plasma is negligible when considering a thermo- 
couple made up of metal and plasma. The 
thermoelectric power of the plasma is esti- 
mated from an approximate solution of the 
Boltzmann equation in the presence of an elec- 
tric field and a temperature gradient. For 
representative plasma conditions covered in 
the experimental work, the analysis gives a 
thermoelectric power of 1.46 x 107° volt/°K. 
For a temperature difference of 2000°K, this 
gives a Seebeck emf of the order of 3 volts; 
this is the order of magnitude of open-circuit 
voltage observed in the various experiments, 
thus giving reasonable agreement with the 
thermocouple theory used. 

An attempt is made in reference 16 to obtain 
an understanding of the characteristics of the 
plasma. Four principal current-limiting fac- 
tors in the operation of the plasma thermo- 
couple are involved. The factors are as fol- 
lows: 

1. The current-emission limit of the emitter 
electrode. 

2. The ohmic resistance of the plasma, due 
primarily to multiple scattering of electrons 
by the plasma ions. 

3. The space-charge region in the vicinity 
of the cathode surface (the main body of the 
plasma is electrically neutral on the average). 

4. The reduced ion density in the vicinity of 
the anode surface where the gas temperature 
is relatively low. 


Work is continuing at LASL to explore ex- 
perimentally the relative importance of these 
current-limiting factors and to determine more 
fully the performance capabilities of the therm- 
ionic cell. 
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REACTOR PHYSICS: 





H,O-MODERATED CRITICAL ASSEMBLIES AND ANALYSES 





The results of a number of critical experiments 
with water-moderated assemblies have been 
reported recently. Three of these, the experi- 
ments on the Martin Power Reactor,’ the 
Yankee critical experiments,’ and the critical 
experiments for the N.S. Savannah’ core, are 
directed toward the investigation of specific 
reactor designs. They treat the usual quanti- 
ties involved in such an investigation: critical 
masses, thermal-flux distributions, control- 
rod effectiveness, reactivity coefficients, the 
appropriate lattice characteristics, and effects 
of water channels. Aside from their usefulness 
in the experimental design of the reactors in 
question, these investigations add to the general 
body of experimental knowledge of H,O-mod- 
erated reactors and are useful for the general 
understanding of such reactors and for the 
testing of theoretical methods. 

The Martin critical experiment investigated 
mock-ups of two full-scale core configurations, 
both about 22 in. in diameter, as well as 
smaller clean assemblies. The reactor is 
highly enriched, using fuel elements of tubular 
shape, 0.375-in. O.D. and 0.315-in. I.D., com- 
posed of a UO,—stainless steel— B,C (burnable 
poison) cermet clad with 0.006-in.-thick stain- 
less steel. The active core length is 23 in. 
Among the measurements of general interest 
are measurements of the effectiveness of (1) 
rare earths as control-rod materials and (2) 
thermal-neutron- and gamma-flux traverses 
through the relatively thin H,O reflector, steel- 
H,O thermal shield, pressure vessel, and sur- 
rounding H,O. 

The Yankee experiments cover assemblies 
which approximate the Yankee fuel-element and 
control-rod lattice but which do not mock up 
the full-size core. The fuel elements are com- 
posed of UO, pellets, 0.30 in. in diameter, 
jacketed in stainless steel. The active length 
is 49.42 in. Three different lattice pitches 
giving water to equivalent uranium volume 
ratios of 2.19, 2.93, and 3.87 were investigated. 
The measurements of lattice characteristics, 


temperature and void coefficients of reactivity, 
water channel effects, and resonance absorp- 
tion in stainless steel are of general interest. 

No comprehensive report has yet been issued 
on the N.S. Savannah experiments, but a paper 
has been presented at a meeting of the Ameri- 
can Nuclear Society. The fuel is uranium oxide 
of 4 per cent enrichment, and jackets of both 
stainless steel and aluminum have been in- 
vestigated, at water/UO, volume ratios of 1.15 
and 1.54. The experiment mocks up the full- 
size N.S. Savannah core. 

Reference 4 describes a more general type 
of experiment which was carried out to study 
the general properties of seed-and-blanket type 
reactors, and to verify the nuclear design 
techniques for such reactors. The experiments 
were carried out in slab geometry. The H,O- 
moderated seed region was a slab, made up of 
highly enriched uranium-zirconium plates, the 
slab being approximately 7 in. thick, 54 in. 
high, and of variable length up to approximately 
40 in. This slab, which was supercritical, was 
provided with cruciform control rods; blanket 
regions composed of water-moderated natural- 
uranium plates could be placed adjacent to it 
on both its large faces. Reactivity measure- 
ments were made and power distributions were 
plotted for various seed-and-blanket combina- 
tions. Temperature coefficients and lattice 
characteristics were also investigated. 

A peculiarity which is perhaps worth noting, 
although its practical effect in power reactors 
is probably small, is what appears to be 
anomalous behavior of the temperature coef- 
ficients of reactivity for all these assemblies 
in the temperature range near room tempera- 
ture. The Martin coefficients show the most 
regular behavior, but the measured reactivity 
losses with increasing temperature over the 
range from 20 to 60°C are about twice as large 
as those calculated theoretically. In the Yankee 
case, when the coefficients are plotted as 
functions of temperature, the curves appear to 
go through minima, or possibly plateaus, near 


15 
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50 or 60°C, before continuing to decrease with 
higher temperatures. The N.S. Savannah ex- 
periments show positive temperature coeffi- 
cients of reactivity up to temperatures of about 
60°C, and the seed-and-blanket experiments 
show a positive coefficient for the seed slab 
alone, up to a temperature of about 50°C, but a 
negative coefficient over this temperature range 
when the seed slab has blanket slabs adjacent 
to it. Whether or not these rather unexpected 
behaviors are easily understandable in terms 
of the peculiarities of the assembly is not ap- 
parent, but it is possible that some general 
theoretical investigation of the temperature 
coefficient behavior in the temperature range 
below boiling would be worthwhile. 

Reference 5 describes a set of critical ex- 
periments and an accompanying analysis which 
were undertaken for the purpose of making a 
rather rigorous test of the accuracy with which 
present calculation methods can determine 
thermal-neutron fluxes in water-moderated 
assemblies. The assembly consisted ofanearly 
homogeneous slab core completely surrounded 
by water reflector. The core, composed of 
water, zirconium, plastic, and highly enriched 
U;03,, had dimensions roughly 6 by 40 by 31 in. 
Flux traverses—or more accurately, activa- 
tion traverses—were made with bare and 
cadmium-covered detectors along the line which 
coincided with the short axis of thecore. Three 
configurations were investigated: the clean core 
surrounded by reflector; the core with a 1-in. 
water gap bisecting its short dimension, and 
surrounded by H,O reflector; and the core— 
water gap configuration with 0.03-in. cadmium 
plates lining the boundaries of the water gap, 
and the core again surrounded by water re- 
flector. The uranium density in the core was 
relatively high, and consequently thermal-flux 
peaks, of height two to three times the average 
thermal flux in the core, existed inthe reflector 
and in the unpoisoned water gap. The arrange- 
ment thus represented a rather difficult test 
for theory because of the large differences in 
properties of the adjacent regions and the small 
sizes of the regions. 

Since the spatial activation distribution of 
foils is proportional to the spatial thermal- 
neutron-flux distribution only if the neutron 
spectrum is constant over the region investi- 
gated, the comparison between experiment and 
theory was made, not on the basis of neutron 
flux, but on the basis of the spatial distribution 


of foil activation. The neutron flux was calcu- 
lated theoretically, and for each foil type 
(manganese, gold, U***) the subcadmium activa- 
tion distribution was calculated, using the theo- 
retical thermal-neutron spectrum, for com- 
parison with the experimental distributions. 

The following calculational techniques (plus 
some corrections for the small degree of in- 
homogeneity which existed in the core region) 
were the major ones involved in the theoretical 
determination of the activation distribution: 

1. A four-group diffusion theory calculation 
of the epithermal neutron-flux distributions 
was used in order to arrive at the spatial dis- 
tribution of thermal-neutron source strength. 
For this calculation the epithermal constants 
of the assembly regions were calculated by an 
infinite-medium, multigroup P-1 code (MUFT). 
Since there was some question as to the validity 
of the infinite-medium constants in the narrow 
regions involved, check calculations of the 
thermal-neutron source distribution were made 
using a P-1 multigroup code which did not 
assume separability of space and energy. The 
difference between the two methods amounted 
to only 1 per cent in the calculated distribu- 
tions. 

2. The thermal-neutron constants for the 
several regions were taken as the Wigner- 
Wilkins averages over infinite mediums of the 
same composition (SOFOCATE code). 

3. The transport of thermal neutrons was 
approximated both by the double P-1 approxi- 
mation and by the diffusion approximation. 

4. Once the thermal-neutron-flux distribu- 
tion was determined, a problem remained in 
determining the activation distribution because 
of spectrum differences from region to region. 
For example, the subcadmium activation of a 
1/v absorber is proportional to the total num- 
ber of neutrons per cubic centimeter in the 
subcadmium energy region, rather than to the 
subcadmium neutron flux. Far from the core- 
water interfaces, the average neutron velocity 
(v), as determined by the Wigner-Wilkins cal- 
culation, was 20 per cent greater in the core 
region than in the reflector region. To deter- 
mine av for the calculation of activation in the 
region near the interface,? was assumed to 
vary linearly between the two values, over a 
distance of one diffusion length on each side of 
the core-reflector boundary. 

The agreement between measured and cal- 
culated activation distributions was gratifying 
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and leaves little doubt that the calculational 
methods tested are entirely adequate for re- 
actor design calculations involving situations 
comparable to those investigated. Even the 
diffusion-theory results are quite good interms 
of the accuracies needed for design. It appears 
that only in special cases, or in connection 
with fundamental neutron-transport studies, 
would further refinements on the methods— 
such as a point-by-point variation of the neu- 
tron spectrum—be called for. 

Reference 6 reports an investigation of the 
neutron-spectrum characteristics of a small 
subregion embedded in a reactor core of dif- 
ferent composition, and particularly of the 
degree to which measurements of the micro- 
scopic core constants in the subregion can be 
considered to duplicate those of a core having 
a composition the same as that of the sub- 
region. The conclusion is that such measure- 
ments are possible and valid and that the 
scheme offers a means of investigating micro- 
scopic core parameters over a wide range of 
core compositions without the use of large 
quantities of experimental fuel elements. 

The subregions investigated had radii of 5.8 
and 9.0 cm and were located at the center of a 
core having a total radius of approximately 
30 cm. Both core and subregion were water 
moderated; slightly enriched uranium and UO, 
were used as fuel alternatively for the core 
and/or the subregion. 

Theoretical calculations were made by a 
four-group method and indicated that, for the 
larger subregion, the ratios of flux densities 
of the four energy groups maintained a nearly 
constant value out to a radius of 5 or 6 cm. 
This indicates, of course, that the neutron- 
energy spectrum is nearly constant near the 
center of the subregion and approaches very 
closely the equilibrium distribution which would 
characterize a complete core having the com- 





position of the subregion. These are the only 
requirements for the validity of the method. 
In the experimental part of the investigation, 
fast and thermal activation distributions were 
measured through the subregion and the sur- 
rounding core and were found to check rather 
closely with the calculated distributions. Ther- 
mal utilization, resonance escape probability, 
and the fast-fission effect were measured in 
the subregions of two different sizes and also 
in a core having the same composition as the 
subregion. The subregion results checked very 
closely with the full-core results for all quan- 
tities except resonance escape probability, 
where there was a discrepancy of 3 or 4 per 
cent in (1—~). This was thought to be con- 
sistent with the small spectrum variations 
indicated by the theoretical calculations. 
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Collapse of Parallel-plate 
Fuel Assemblies 


The design of high specific power reactors in- 
volves choosing a coolant mass flow rate thatis 
as high as possible, consistent with an acceptable 
pressure drop and with a continuing structural 
integrity of the fuel element. With respect to 
this latter item, reference 1 presents an ana- 
lytical study of the flow velocities necessary to 
collapse parallel-plate fuel assemblies, suchas 
are found in MTR-ETR fuel elements. The 
assumptions used in the derivations are as 
follows: 


1. The plates are homogeneous, isotropic, elastic, 
initially flat or uniformly curved, uniform in spacing 
and dimensions, and free of unidentified sources of 
deformation. 

2. The coolant is incompressible; all channels 
have the same mass flow; at any cross section 
normal to the longitudinal axis, the flow within any 
channel is uniform; and leakage between channels 
is suppressed. 

3. Plates are broad enough (in comparison with 
their thickness) so that shear deformation is negli- 
gible, and they are long enough (in comparison with 
their breadth) to enable them to deflect locally 
without significant redistribution of flow among the 
coolant channels, 

4. Side plates or supports are rigid. 


The various cases and resulting formulas are 
presented below: 
Case A: Flat Plates with Built-in Edges 


l, 
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Case B: Flat Plates with Simply Supported 
Edges 
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Case C: Flat Plates with Longitudinally Spaced 
Supports 
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Case D: Curved Plates with Fixed Edges 
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Case E: Curved Plates with Hinged Edges 
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where Vc =linear velocity of coolant for 
‘which collapse occurs due to 
flow-induced effects alone, ft/sec 
g = gravitational acceleration, ft/sec” 
E = Young’s modulus, lb/sq ft 
a = plate thickness, ft 
h = initial flow channel thickness at 
midspan, ft 
p = coolant density, lb/cu ft 
b = width of flat plate or chord of 
curved plate, ft 
v = Poisson’s ratio 
1 = longitudinal distance between 
plate-edge supports, ft 


2 
By =k fe (a) +f; (a) 


B. =—— fz (a) + fy (a) 


R = initial mean radius of curvature, 
ft 
IJ =moment of inertia of beam cross 
section per unit width of beam, 
ft‘/ft 
A =area of beam cross section per 
unit width of beam, ft?/ft 


1 a sin’ a 
Ix (@) = > Sin? a = a 
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fa (a) = f,(0) — 35 


1  — a@ cos’ q 
fs (@) = 5 Gin? a g (2 ~ sin 2a + 2S a ) 


2a =curved-plate arc between supports 


The effect of different pressures on either 
side of the plate is discussed, as is the case 
where the fuel plates are subjected to external 
transverse loadings in addition to hydraulic 
forces. Some experimental justification for the 
approach is presented in reference 1, but the 
author recommends additional experiments to 
validate the analytical study. 

Additional data on collapse of parallel-plate 
assemblies are presented in reference 2. Plexi- 
glas parallel plates were tested in a flow 
facility wherein the side plates were rigid. 
The Case A equation predicted a collapse 
velocity of 21 ft/sec. Results are as follows: 


Flow through the test section was gradually in- 
creased until plate-to-plate contact was seen. Slight 
plate contact was visible at 21 ft/sec; full contact 
occurred at 31 ft/sec at the leading edge of the 
plates. At 41 ft/sec, contact extended downstream 
and gradually tapered off until it ceased 4 to 5 in. 
above the inlet edge of the plates. A definite flow 
disturbance existed across the plates and was made 
visible when air was injected into the flow stream. 


The work of reference 1 is also extended to 
include the effect of flow redistribution on the 
critical flow velocity. 


Nucleate Boiling 


In the design study’ of the large organic- 
cooled and -moderated reactor which was made 
as part of the AEC survey of large power reac- 
tor types under Public Law 85-590 (reviewed 
in Power Reactor Technology, 2(4): (September 
1959)), one of the anticipated technological ad- 
vances was the use of nucleate boiling of the 
organic coolant in the hot-coolant channels. 
Since the thermal performance of the design 
depended rather strongly on this feature, it 
would appear that the question of nucleate 
boiling in organic coolants is one of consider- 
able importance. Reference 4 concerns nucleate 
boiling characteristics under conditions of pool 
boiling at atmospheric pressure. The following 
coolants were investigated: 





1. Diphenyl. 

2. Metaterphenyl (Santowax-M). 

3. Orthoterphenyl (Santowax-O). 

4. Mixture of 4.35 wt.% paraterphenyl in 
metaterphenyl. 

5. Eutectic mixture of 32 wt.% diphenyl in 
orthoterphenyl. 


The data were presented in graphical form 
and will not be reproduced here. A peak heat 
flux was found for all coolants studied except 
orthoterphenyl. 

The case of pool boiling is not directly appli- 
cable to organic-cooled power reactors since the 
reactor coolant would be under conditions of 
forced convection, probably at elevated pres- 
sure, and subject to varying degrees of sub- 
cooling. Hence the authors presented methods 
of estimating the effects of pressure, velocity, 
and subcooling as well as compared their data 
with other pool boiling investigations. 


The effects of velocity and subcooling were 
examined by using the equations of Core and 
Sato> for diphenyl and Santowax-R extrapolated 
to zero velocity. The results of the present 
study showed peak fluxes of 95,000 Btu/(hr) 
(sq ft) for diphenyl and from 97,000 to 135,000 
Btu/(hr)(sq ft) for the other substances; extrap- 
olations of Core and Sato’s equations yield 
values of 116,000 Btu/(hr)(sq ft) for diphenyl 
and 152,000 Btu/(hr)(sq ft) for Santowax-R. The 
lowest diphenyl velocity used by Core and Sato 
in their experiments was 0.5 ft/sec, whereas 
that for Santowax-R was 5 ft/sec. 

The correlations of Rohsenow® and Forster 
and Greif’ for pool boiling of other liquids were 
examined as to their applicability for use with 
diphenyl, orthoterphenyl, and metaterphenyl. 
The Rohsenow correlation was found to be 
satisfactory for diphenyl and metaterphenyl, but 
the correlation did not hold for the case of 
orthoterphenyl. One of the correlations sug- 
gested by Forster and Greif was found to be 
satisfactory for diphenyl and metaterphenyl, but 
it too was unsatisfactory for use with orthoter- 
phenyl. Several coefficients used inthe correla- 
tions were determined. 

An interesting phenomenon was observed 
during the runs with diphenyl. During one of the 
runs, the test section was allowed to film boil 
for 15 to 30 sec, and a thin film of amorphous 
carbon deposited on the heater. After the test 
section was brought back to the nucleate boiling 
regime, the heat-transfer coefficient showed 
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marked improvement. Figure 8 shows the im- 
provement after film boiling. These data were 
found to be reproducible in test runs conducted 
over a year apart. It was concluded by the 
authors that the thin layer of carbon so improved 
the nucleation mechanism that, for a given heat 
flux, the temperature difference requirements 
were lowered. A pertinent question concerns the 
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Figure 8—Comparison of heat flux versus tempera- 
ture difference for diphenyl, before and after film 
boiling.‘ A, diphenyl after boiling. B, diphenyl before 
boiling. 


significance of the observations to organic- 
cooled power reactors employing boiling. One 
point is that the modified surface appears to 
be more stable to nucleate boiling, since it 
appears that the test section did not film boil 
at heat fluxes well over the peak heat flux 
recorded for the clean surface. Thus an organic- 
cooled reactor design that is power limited by 
the onset of film boiling might profit if element 
surfaces could be modified to give the heat- 
transfer coefficients implied by the upper curve 


of Fig. 8. A further question, of course, is 
whether the carbon deposit, helpful in a thin 
layer, might be harmful in thicker layers, and 
whether the presence of nucleate boiling might 
favor the build-up of deposits over along period 
of time. 

The last issue of this publication (Vol. 2, 
No. 4, September 1959) discussed the role of 
surface condition in boiling water and cited 
Argonne National Laboratory (ANL) data which 
revealed that surface irregularities tend to 
increase the peak heat flux because of increased 
turbulence and more numerous nucleation points. 
References 8 and 9 are reports of recent re- 
search in the area. The reports will not be 
reviewed in detail since they are of primary 
interest to the specialist in boiling heat transfer. 
Reference 9 describes a new approach to the 
determining of the number of active sites present 
during nucleate boiling by simultaneously elec- 
troplating a layer of metal on the boiling surface 
during boiling and subsequently by counting the 
pinholes where no metal was deposited. The 
conclusions of the two reports are believed to be 
of sufficient general interest that they are 
reproduced herein; those of reference 8 are as 
follows: 


1. Nucleation occurs from pre-existing gas-filled 
cavities on the surface. 

2. A single dimension in general is sufficient to 
characterize such an active cavity. 

3. The wall superheat at which such a cavity will 
become active is fixed by surface tension and other 
known fluid properties. 

4. The mean surface temperature is not equal to 
the temperature felt by an active cavity. 

5. For surfaces made of the same material and 
treated in the same way,a single plot of the number 
of cavities versus a group with the dimensions of 
length can be developed which is appropriate to 
apply to a variety of fluids at various pressures. 


The conclusions of reference 9 are: 


1. The relation between heat flux and active site 
population is 


Q/A = 1400 (N/A)? 


where Q/4 = heat flux, Btu/(hr)(sq ft) 

N/A = number of active sites on heated sur- 
face per square foot of heat-transfer 
area 

2. At active site populations greater than 500 per 
square foot, the heat-transfer coefficient, h, can be 
represented by 
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h = 49 (N/A)?*48 


3. The amount of heat attributed to each active 
site becomes less as the population of active sites 
increases, At the beginning of boiling, although 
most of the heat is transferred by free convection, 
the first active site in the 1.77-sq in. test area (or 
81 sites per square foot), was associated with 
11,000 Btu/hr. By comparison, at a heat flux of 
317,000 Btu/(hr)(sq ft), each active site was asso- 
ciated with only 1.95 Btu/hr. 

4, Electroplating during boiling is a new tool for 
boiling research. The populations of active sites on 
a heated surface can now be determined at much 
greater heat fluxes than possible by previous tech- 
niques. An analysis of information from additional 
experiments could lead to a unified theory for the 
entire range of nucleate boiling in which the surface 
condition, as well as nucleation and bubble growth, 
is treated as a variable. 


Liquid-metal Heat Transfer 


Reference 10 is an analytical study of turbulent 
liquid-metal heat transfer in channels. The case 
considered is restricted to the entrance region 
of a channel, where boundary layers are not, by 
definition, fully established. The idealized sys- 
tem investigated is the following: 


1. The wall temperature distribution is a simple 
step function, ! =t) for x <0 and ¢=¢t, for x >0. 
The fluid approaching the entrance region has an 
established temperature, ¢p. 

2. Longitudinal heat conduction is small com- 
pared to convection and is neglected. 

3. The established turbulent velocity profile is 
represented by a uniform distribution uw =U. 

4. The eddy diffusivity distribution varies linearly 
with distance from the wall and as the nine-tenths 
power of the Reynolds modulus, ¢€/v = C,Re’*(y/b). 

5. The fluid properties are invariant with tem- 
perature. 

6. Steady state exists. 


€ = eddy diffusivity, sq ft/hr 
v = viscosity, sq ft/hr 
C; = dimensionless constant 
Re = Reynolds number 
b = channel half-thickness, ft 
y = distance measured into channel from 
wall perpendicular to the channel 
boundary, ft 


Here 


Several of the above postulates bear discussion. 
The effects of axial conduction anda nonuniform 
velocity profile were investigated and found to be 





unimportant over the range of variables studied. 
The diffusivity function shown was formulated 
under the assumption that the heat-transfer 
and momentum-transfer eddy diffusivities are 
identical, and the postulated function is a linear 
approximation to the actual momentum-transfer 
eddy diffusivity. Since the actual momentum- 
transfer eddy diffusivity is constant in the 
turbulent layer, the postulated function departs 
radically from the actual function as distance 
into the turbulent layer is increased; however, 
the postulated function is a good approximation 
in the laminar and buffer layers. It was calcu- 
lated that the Nusselt moduli obtained with the 
linear function were within 10 per cent of the 
moduli when the actual momentum-transfer eddy 
diffusivities were employed in the calculations. 
The linear approximation makes the equations 
more tractable, and the resulting answers are 
of sufficient accuracy because the heat-transfer 
limitations are not set by the turbulent layer. 
The differential equation for the convective 
transfer of heat was solved by a separation-of- 
variables technique, and the symmetry boundary 
condition resulted in the typical infinite series 
solution, since the differential equation yielded 
zero-order Bessel functions of the first and 
second kind for its general solution. The eigen- 
values and series coefficients were evaluated 
for a range of Reynolds and Prandtl moduli, 
and Nusselt numbers and mixed-mean fluid 
temperatures were calculated. Prandtl numbers 
from 0.002 to 0.06 and Reynolds numbers from 
5000 to 1 x 10° were used. Figure 9, typical of 
the results, illustrates the fact that entrance 
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Figure 9—Nusselt modulus versus X and Re;Pr = 
0.002 (reference 10). 
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effects are present for many channel thicknesses 
downstream, particularly at high Reynolds num- 
bers. The calculated results were compared with 
other work, both experimental and analytical, 
and good agreement was obtained. 


Liquid Metals with Internal ° 
Heat Generation 


The system studied in reference 11 consists 
of a long, smooth, insulated pipe with liquid 
metal flowing under internal heat-generation 
conditions. Guard heaters were used around the 
insulation to achieve adiabatic operation. 

A general theory of heat transfer with internal 
heat generation was derived, and experiments 
were conducted to validate the theory. The 
solution for the conservation-of-energy equation 
yielded an analytical expression for the following 
parameter T as a function of Np, and Np;: 


k (Tw — Tm) =T 
Im Ww 
where & = fluid thermal conductivity, Btu/(hr) 
(ft)(° F) 
T, = wall temperature at fluidinterface, 
°F 


Tm = mixed mean fluid temperature, °F 

9m = mean volumetric heat generation 
rate, Btu/(hr)(cu ft) 

Yw = pipe radius, ft 


The analytical expression was derived under 
the following assumptions: 


1. A steady state exists. 

2. Temperature and velocity distributions are 
independent of axial position and symmetrical about 
the axis. 

3. Constant radial heat flux at the wall. 

4. Heat sources exist in the fluid and are inde- 
pendent of axial position and symmetrical about the 
axis. 

5. Physical properties are not functions of tem- 
perature, 

6. Physical properties are not affected by con- 
ditions of fluid flow. 

7. An analogy exists betweenheat and momentum 
transfer for turbulent flow; i.e., €, = @€m. 

8. Heat is transferred by convection, molecular 
conduction, and eddy conduction. 


Since one of the purposes of the paper was to 
study the relations between €, and €,,,, the deri- 
vation employed a value of unity for @. 


The experiments were conducted with mercury 
for Npr from 0.0191 to 0.0224 and Nre from . 
29,000 to 164,000. The fluid was internally 
heated electricaily, and values of T were de- 
termined by measuring appropriate tempera- 
tures and other necessary parameters. 

The experimental values of T, Texpt were 
found to be on the average 49 per cent higher 
than the values calculated using the a = 1 anal- 
ogy. Based on this result, the author concludes: 
“Thus, it may be concluded that, onthe average, 
the radial conductances are smaller than those 
calculated from the heat-momentum-transfer 
analogy. This is essentially the same result as 
ordinary heat-transfer tests with liquid metals 
have shown. Thus it appears that the ratio a is 
not unity but somewhat smaller.” 

The difficulty with the experiment lies in the 
fact that the modulus T is calculated by meas- 
uring the difference between two large numbers 
of nearly the same value, 7, and T,,. Typical 
differences of 2 to 6°F were measured, and 
values of wall temperature ranged from about 
100 to nearly 200°F. This difficulty is reflected 
by the fact that the value of Texp¢ showeda 
maximum of +35 per cent variation from the 
mean. Temperature differences were stated to 
be accurate to +40 per cent at the worst. 


Water Flowing Parallel 
to Tube Bundles 


The fuel element consisting of bundles of fuel 
rods cooled by water which flow parallel to the 
axes of the rods is of continuing interest to the 
reactor design engineer. Recent experimental 
and analytical data’? have shown that the Colburn 
equation may not be correct for determining 
turbulent, nonboiling heat-transfer coefficients. 

The recommended equation is as follows: 

Ae =C (Re)’* (Pr)' 
where h =heat-transfer coefficient 
De = 4 X hydraulic radius 
k =thermal conductivity based on film 
temperature 
Re = Reynolds number based on film 
temperature 
Pr = Prandtl number based on film tem- 
perature except Cy , which is based 
on bulk temperature 
C = empirical constant 
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Values of C are given by the following equations 
for 2.5 x 10‘ < Re <10°: 
Triangular pitch, 


C = 0.026 (S/D) — 0.006 for 1.1< S/D < 1.5 


where S = center-to-center distance between 
tubes 
D =tube diameter 


Square pitch, 
C= 0.042 (S/D) — 0.024 for 1.1 < S/D < 1.3 


This new correlation usually yields higher heat- 
transfer coefficients than does the Colburn equa- 
tion; at large values of S/D, the improvement 
is about 40 per cent. Also, the square lattice, at 
a given S/D, shows improved coefficients com- 
pared to the triangular lattice. 
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Delayed Neutrons 
in Reactor Kinetics 


Although extreme accuracy in the calculation of 
kinetic behavior is not often required for reac- 
tor design purposes, many reactor experiments 
depend, for the measurement of reactivity, on 
kinetic measurements. The most common of 
these is the measurement of an asymptotic 
positive reactor period in the range from a few 
seconds to a few hundred seconds. In such cases 
an accurate knowledge of the kinetic character- 
istics of the reactor is essential for accurate 
results. Three main types of effects must be 
taken into account: uncertainties in the delayed- 
neutron characteristics for a given fissionable 
isotope, the effect of neutron importance in 
determining the effective values of the delayed- 
neutron fractions, and delayed-neutron produc- 
tion by sources other than the main fissionable 
isotope (e.g., by fast fission in U?** or by y-n 
reactions). Reference 1 treats the first two of 
these questions. 

The value of a delayed neutron to the chain 
reaction is, in general, greater than that of an 
average prompt neutron because the delayed 
neutrons are produced at lower average energies 
than the prompt neutrons and, therefore, have a 
lower probability of leaking from the reactor. 
For this reason the effective delayed-neutron 
fractions may be larger than the true fractions. 
This effect is small in large low-leakage reac- 
tors but may be appreciable for small reactors, 
particularly those moderated by H,O, in which 
the slowing-down length is quite sensitive to 
neutron birth energy. The kinetic effects of 
delayed-neutron importance can be discussed 
conveniently in terms of the inhour equation, 
which relates the asymptotic reactor period to 
the reactivity. 

For long periods the inhour equation 


1* B; 
ee ee 
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can be reduced to 
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where p is the reactivity, /* is the effective 
prompt-neutron lifetime, 7 is the asymptotic 
reactor period, A; is the decay constant of the 
ith species of delayed-neutron precursor, and 
B; is the effective fraction of total fission neu- 
trons resulting from the 7th precursor species. 
The £; is, in general, different from 8;, the true 
ith delayed fraction, but is, of course, dependent 
on £;. 
The reference rewrites Eq. 1 in the form 


-_ 4/8 
p = lB len » I; (7) (2) 


where £ is the true total delayed fraction and 
I’;(T) is given by 


B,/ 
aM eat 
Thus the two terms, 8 and XI; (tT), depend only 
on the true delayed-neutron fractions, and all the 
effects stemming from the importance of the 
delayed neutrons in the reactor are lumped into 
the term [8/8] (7), defined as 


7 d Ete, (7) 
B Ion _ 7 Bi 

B 7 Tr; (T) 

It is thus an average of the (effective 6;)/(true B;) 
weighted by I(T). It is evidently a function of T, 
as indicated, and as is to be expected, since the 
relative importance of the various delayed 
groups in determining 7 varies as 7 varies. The 
term [8/8](T) would be expected to be insensitive 
to variations in the experimental values of 8;, 
although sensitive to the energies of the delayed 
groups. Therefore Eq. 2 has separated the terms 
which (1) involve experimental uncertainties in 
8B; and A; and (2) are common to all reactors 
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employing a given fissionable isotope, from the 
term which (1) involves the delayed-neutron 
importance, (2) is subject to experimental error 
only with respect to delayed-neutron energies, 
and (3) will vary from reactor to reactor. 

Reference 1 investigates the effective delayed- 
neutron fractions for a case in which the ratios 
B; /B; would be expected to be relatively large: 
a thin (7-in.) highly enriched water-moderated 
slab, reflected by water, and the same slab, re- 
flected by a natural-uranium and water blanket. 
The ratio B;/8, was determined for each de- 
layed group by making four-group criticality 
calculations in which the effects of the 7th group, 
with its appropriate birth energy, was treatedas 
a separate fast-neutron source. From the ratios 
so obtained, the quantity [8 /8] (rT) was computed 
for various values of tT from 20 to 300 sec. The 
values obtained for the H,O-reflected slab were 
all approximately 1.26, and the variation over 
the range of 7 was no greater than the uncer- 
tainties in birth energy of the delayed groups. 
It may therefore be concluded that an effective 
value of the total delayed fraction, 8, can be 
derived for a given reactor and that the ratio 
B/B can be used as a correction term for the 
inhour equation using true values of §;, to take 
account of delayed-neutron importance over the 
range of periods normally used in experiments. 

The reference also tabulates, for various 
values of 7 from —5000 sec to +200 sec, the 
values of TI; and p = I;as derived from three 
different sets of experimental data: the Hughes 
data’ of 1948, the Keepin data® of 1955, and the 
Keepin data‘ of 1957. 

Reference 2, which is a rigorous general dis- 
cussion of the application of reactor kinetics as 
a tool in experimental reactor physics, is also of 
interest in connection with this subject. 


Reduced Delayed-neutron 


Group Representations 


Reference 5 reports an analytical study ofthe 
validity of representing the six delayed-neutron 
groups by a smaller number of groups in kinetic 
calculations. The kinetic behavior as calculated 
with “equivalent” one-group, two-group, and 
three-group representations are compared with 
those calculated using the six actual delayed 
groups for a variety of different cases. These 
cases include the behavior of the asymptotic 


reactor period (as given by the inhour equation), 
the reactor transfer function, and the behavior 
after step increases and step decreases of 
reactivity. It is shown that the two-group and 
three-group approximations duplicate reason- 
ably well the behavior with the full six groups. 
The principal advantage of the three-group 
representation over the two-group is that, by 
choosing the characteristics of the third group 
to duplicate those of the actual delayed group 
having the longest half life, the asymptotic rate 
of decay of power in the subcritical reactor can 
be duplicated. 


Approximations for Reactor 


Accident Calculations 


Reference 6 discusses methods of approxi- 
mating the behavior of reactor power during a 
typical start-up accident or during similar 
transients. The typical assumed accident is one 
in which reactivity is added continuously to the 
initially subcritical reactor (usually at a constant 
rate) until the reactor is made subcritical by the 
over-power scram circuit or by other means. 
In the reference the accident is divided into 
three stages: the increase of reactivity to de- 
layed critical, the interval from delayed critical 
until the shutdown mechanism begins to function, 
and the shutdown period. Different approxima- 
tions, which are not amenable to brief summary, 
are applied in the three stages. 

Although accidents of this type canbe investi- 
gated to any desired degree of numerical accu- 
racy by machine computation, the uncertainties 
in the physical processes involved do not often 
justify very great accuracy of computation. 
Consequently, the approximations discussed in 
the reference may be found to be useful for 
many practical cases as well as helpful for the 
insight they give into the kinetic behavior of 
the reactor. 

Reference 7 also treats approximate methods 
for calculating the effects of hypothetical reac- 
tivity accidents. Specifically, it covers the cal- 
culation of rates of reactivity removal necessary 
to limit reactor power to a specified value when 
the power is rising on a period T because of the 
step insertion of an excess reactivity 5k». 

For dkp < 8, with reactivity reduced at acon- 
stant rate y sec’ after time ¢ = 0 and power 
Py, the power P at any time ¢ > 0 is given by 








26 GENERAL RESEARCH AND DEVELOPMENT 


* 
cae wr |! 52d.) ery 
Py B 


272 
3 (1 — alee r ‘ -(* = ca _ e-Bi/I*) 





2] * 2 
4 aA i. Fe bkp 42 
B 2p 


where @ =y/(Adkp), l* is the effective prompt- 
neutron lifetime, 8 is the delayed-neutron frac- 
tion, and A is the decay constant of a single 
group of delayed neutrons equivalent to the six 
actual delayed groups. The reference states 
that an approximation which will give conserva- 
tive results (i.e., a high value of P) isA= 
= A,B;/B. The equation above is based on the 
linearized kinetic equations; consequently, it 
should be valid only if y is chosen large enough 
that the maximum value of P is not greatly 
larger than P). It is stated that values of y which 
will reduce the reactivity to zero in about 2% sec 
or less will limit P to values not substantially 
greater than Py, provided 5k» < 8B. 

In the case of 5kp > Bthe power will rise 
above the value P, unless the reactivity can be 
reduced nearly instantaneously. If reactivity is 
reduced at the rate y sec™!, the maximum power 


P,, is approximated by 
Pm _ (6kp — B)’ 
lo pM LR aint fees At Ma 
Be Py 2yl* 


If reactivity is reduced at an accelerated rate a, 
where a has the units (5% of¢/Reps )/sec”, the maxi- 
mum power is given by 


rm Pm _ 2 V2 (dkp — p)” 
be Py 31 va 





The reference compares the results of the 
above approximations with digital computer re- 
sults, using the exact space-independent kinetics 
equations to demonstrate their validity. 


Oscillation Experiments on 
Large Heavy-water Reactors 


Oscillation experiments and their analyses to 
determine temperature coefficients of reactivity 
for a large, natural-uranium reactor, cooled and 
moderated by D,O, are reported in reference 8. 
The experiments are interesting primarily for 


the techniques they establish, especially since 
they were made on a large reactor whichhad no 
special provisions for reactivity oscillation. 

A trapezoidal reactivity oscillation was used, 
produced by periodically withdrawing a gang of 
control rods at constant speed for a time 
interval ¢ (the drive time), holding the rod posi- 
tion constant for a time interval ¢, (the wait 
time), driving the rods in again at the same 
constant speed for a time interval ¢,, and again 
holding position constant for a time interval fp. 
Values of 2 and 4 sec were used for ?,;, and 
values of 1, 2, 4, 6, and 8 sec were used for /,. 
The reactivity coefficients were derived from 
the observed behavior of the neutron-flux level 
during the latter three-fourths of the wait time. 
The observed behavior was taken as the average 
for 10 to 15 consecutive oscillations. 

The purposes of the wait time and of the 
variations in wait time are to separate the 
effects on the neutron detector of changes in 
flux shape resulting from the displacement of 
the control rods from the effects onthe detector 
of changes in over-all power resulting from the 
reactivity transients brought about by changes in 
the temperatures and in the control-rod posi- 
tions. Such a separation is necessary in a large 
reactor. 

The temperature coefficient was derived by 
fitting the observed flux-level variation with a 
mathematical model in which the temperature 
coefficient of reactivity was treated as an ad- 
justable parameter. In the particular reactor 
investigated, the moderator temperature coeffi- 
cient of reactivity did not appear in the frequency 
range investigated because of the long time 
constant associated with the coefficient. Ap- 
parently, the coolant coefficient of reactivity 
was small; the observed coefficient, which was 
negative, with a delay of 2.3 sec, was attributed 
to the fuel temperature. 


Transfer Functions of 
Distributed Parameter Systems 


In reference 9, transfer functions are derived 
for a solid-fuel reactor—as approximated by 
two different models—and for a counterflow 
heat exchanger. The reactor and heat exchanger 
are considered as distributed-parameter (rather 
than lumped) systems in so far as heat transfer, 
flow, and temperature distributions are con- 
cerned, although the power distribution in the 
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reactor is considered to be constant. Lumped 
network approximations to the distributed- 
parameter transfer functions are given in an 
appendix. 
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AEC Reactor Safety Program 


The U.S. Atomic Energy Commission has is- 
sued a document (reference 1) which outlines 
the comprehensive research and development 
program of the AEC on reactor safety and con- 
tainment. The document also summarizes 
briefly the results that have been obtained in 
the several projects which make up the pro- 


gram and illustrates these results with some 
of the data that have been obtained. Since it is 
important to the nuclear plant designer to fac- 
tor into his reactor and containment designs 
the most up-to-date safety and containment 
information, Table VI-1 has been prepared, 
summarizing the program, indicating the agen- 
cies carrying out the various projects, and 
repeating most of the references from the 
AEC document. 


Table VI-1 AEC REACTOR SAFETY PROGRAM 





Selected 
references 
(from ref. 1) 


Program Main activity or objective Status Agency 





Reactor experiments on Dynamic measurements on experimental 


kinetics and control reactors , 
SPERT-I Atmospheric pressure; natural-circula- Operating Phillips 2-5 
tion H,O; highly enriched core since Petroleum 
1955 Co. (NRTS) 
SPERT-II Pressure to 300 psi; H,O or D,O Critical Phillips 6,7 
with circulation to 20,000 gal/min; 1959 Petroleum 
various cores and reflectors Co. (NRTS) 
SPERT-III Pressure to 2500 psi; circulating Critical Phillips 6,7 
H,O for power up to 60 Mwi(t); December Petroleum 
various cores; boiling or sub- 1958 Co. (NRTS) 
cooled operation 
KEWB Homogeneous water boiler (H,O); Operating Atomics 8,9 
graphite reflected; normal since Inter- 
operating power 50 kwit) 1956 national 
EBR-I NaK-cooled fast breeder; normal Operating Argonne 10-14 
operating power 1.4 Mwit), since National 
250 kw(e); reactor dynamics with 1951 Laboratory 
various cores by means of (NRTS) 
oscillator tests 
Safety-system performance Measure practical limits of safety- In Oak Ridge 
system response and develop progress National 
reliable fast-acting systems Laboratory 
which approach these limits; 
work to date primarily for pool 
type reactors 
Development of self- Develop self-actuating reactor In Atomics 15-17 
actuating safety devices shutdown devices for research progress Inter- 
(reactor fuses) and power reactors national 
Theoretical reactor Develop and test analytical pro- In Phillips 
kinetics cedures; identify and understand progress Petroleum 
basic sources of reactivity change; Co. and 


Ramo-Wool- 
dridge Corp. 
(SPERT-I 
analysis); 
Atomics 
International 


advance technique of safety 
calculations 
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Table VI-1 (Continued) 
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Selected 
references 
Program Main activity or objective Status Agency (from ref. 1) 
Energy and fission-product Predict realistically the form and 
release caused by fuel quantity of fission-product release 
melting following core meltdown; obtain 
‘ understanding of factors influencing 
initiation, propagation, and rate of 
metal-water and metal-gas reactions 
Fission-product release Species and quantity of fission products In Oak Ridge 18-21 
upon melting of irradiated released as affected by nature of progress National 
fuel samples fuel element, concentration of Laboratory 
fission products in element, 
atmosphere of meltdown, time at 
melt temperature, agitation and 
disposal of melt 
Meltdown of zirconium- Investigation of effects of melting In MSA 
uranium subassemblies subassemblies of plates in steam progress Research 
atmosphere; fission-product release, Corp. 
chemical reaction, and physical 
distribution of melt and fission 
products 
Field release of fission Release and dispersion of fission First Convair, 
products from melted products in the field from melted tests Fort 
fuel fuel in 1958, Worth 
to be (NRTS) 
reported 
in 1959-60 
Transient reactor test Meltdown of fuel elements and core TREAT Argonne 22 
facility (TREAT) segments in a special transient Critical National 
reactor to study melting behavior, February Laboratory 
fission-product release, and 1959 (NRTS) 
chemical reactions; initial program 
directed at fast-reactor safety 
Metal-water reaction Studies of chemical reactions between In Argonne; 23-27 
studies water and molten metals of types progress General 
used in reactor construction and the Electric 
rates and magnitudes of resulting (Valle- 
energy releases; theoretical, citos); 
laboratory, engineering, and Stanford 
reactor investigations Research 
Institute 
Metal-gas reaction studies Fundamental study of ignition and In Argonne 28 
burning of metals in gaseous progress National 
atmospheres, work presently con- Laboratory 
centrated on uranium and its alloys 
Containment studies Develop methods for determining the 
loads which might be generated in a 
reactor accident and for designing 
containment structures which will 
resist these loads 
Energy sources Present work includes experimental In Armour 29,30 
study of explosive decompression progress Research 
of hot water (ARF), external effects Foundation; 
of pressure vessel rupture (APG), Aberdeen 
and methods of simulating reactor Proving 
energy releases (SRI, NOL) Ground, 
Stanford 
Research 


Institute; 
Naval 
Ordnance 
Laboratory 
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Table VI-1 (Continued) 





Program 


Main activity or objective 


Status 


Agency 


Selected 


references 
(from ref. 1) 





Pressure transmission 
and reflection 


Missiles 


Earthquakes 


Primary containment 


Blast shields 


External containment 
structures 


Study of hydromechanics of blast- 
wave propagation in various media; 
mathematical description of 
pressure-volume characteristics; 
experimental crushing of materials 
under shock loads; velocity of sound 
in gas-water dispersions 

Study of penetration of materials by 
rod type missiles 


Study of possible effects of earth- 
quakes on reactors 


Experimental investigation of the 
mechanics of large deformation 
and rupture of pressure vessels 

Study of use of blast shields to reduce 
loading on outer containment 
structure; experimental investiga- 
tion of dynamic strength of 
crushable materials 

Investigation of response of outer 
containment structures to transient 
internal loads; establish design 
criteria by measuring elastic and 
plastic response and comparing with 


In 
progress 


In 
progress 


Report 
expected 
1959 


In 
progress 


In 
progress 


In 
progress 


Armour 
Research 
Foundation 


Stanford 
Research 
Institute 

Lockheed 
Aircraft 
Corp. and 
Holmes 
and 
Narver, 
Inc. 

Naval 
Ordnance 
Laboratory 

Armour 
Research 
Foundation 


Aberdeen 
Proving 
Ground 


29,31 


32 


33 


34 


35-37 





theory; verify the usefulness of 


scale-model testing 





Metal-Water Reactions 


The possibility of a reaction of the metals 
contained in a reactor with its coolant influ- 
ences many areas of reactor design. Two 
classes of reactor accidents may be conceived 
in which metal-water reactions might be im- 
portant. The first of these would involve the 
sudden melting of the metal in an environment 
of more or less liquid water, as the result, for 
example, of a rapid power excursion. This type 
of accident, which could conceivably lead to an 
explosive chemical energy release, has been 
considered not applicable in many water reac- 
tor cases. The second class is typified by the 
loss-of-coolant accident, in which the melting 
occurs rather slowly, in a steam atmosphere, 
and molten metal may drop into a pool of re- 
sidual water at the bottom of the reactor vessel. 
Reference 38 illustrates a method for deter- 
mining the extent of a metal-water reaction in 


an accident of this type and applies the method 
to calculate the amount of reaction for a par- 
ticular accident occurring to a particular re- 
actor. The method of treating the problem is 
believed to be of sufficient interest to warrant 
discussion. 

The problem can be divided into two sections: 

1. The application of high-temperature zir- 
conium-water reaction data to determine the 
temperature excursion (and extent of reaction) 
of the reactor core when it is being affected by 
a combination of heat produced by the oxida- 
tion process together with the decay heat. 

2. Determination of a cooling curve, and 
extent of oxidation, for the reaction between a 
molten drop of uranium alloy and/or zirconium 
and steam. This assumes the clad has been 
ruptured, releasing molten drops of fuel. The 
basic method used is that of Lustman.*® 


The kinetics of the reaction of Zircaloy and 
water are contained in the following equations: 
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AW =k(T)t”(7) 
k(T) = Ae-Q/RT 


where W = oxygen weight gain, g/cm’ 
k(T) = temperature-dependent rate constant 
t =time, sec 
n(T) = temperature-dependent exponent of 
the time 
A = constant 
Q = activation energy, cal/mole 
R = gas constant, cal/deg-mole 
T = absolute temperature, °K 


These functions are plotted in Figs. 10 and 11. 
It should be noted that some extrapolation has 
been done in Fig. 10, and this carries over to 
Fig. 11 also. 

The procedure starts with the assumption of 
some sort of credible accident, and this obvi- 
ously is a function of the detailed reactor de- 
sign; the loss-of-coolant accident studied in 
reference 38 is probably one of the most severe 
accidents that can happen to a power reactor, 
since it assumes a rupture somewhere in the 
primary coolant loops such that the core is 
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Figure 11—Metal-Water reaction curves for 


Zircaloy-2.** 


After loss of coolant the temperature history 





left bare with a pool of water standing in the 
bottom of the pressure vessel. Thus water is 
available for a metal-water reaction but not 


for cooling. 
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Figure 10—A(T) and n(T) versus T. Dotted lines are 
extrapolations of data.*® 





is divided up into increments of temperature, 
AT. For each AT the sum of decay and reac- 
tion heats is equated to the heat required to 
raise the core temperature AT degrees, and 
the kinetics at the temperature of the upper 
limit of AT are used to compute the amount of 
reaction over the whole AT range. The basic 
chemical reaction is 


ZY solid + 2H, Ogas — ZrO» solid + 2H gas (1) 


Because Eq. 1 is exothermic, the metal-water 
reaction can. become autocatalytic, heating the 
metal to its melting point (about 2200°K). The 
amount of heat produced by Eq. 1 is 


AHp =—149,120 + 8.57 — 1.44 x 1077 


, 2-41 x 10° 
T 


where AHp= calories released from Eq. 1 as 
written* 
T = temperature, °K 


All the heat released by Eq. 1 does not go into 
the metal, however, since some heat must go 
into heating 2 moles of steam to the tempera- 


(2) 





* Reference 38 shows the units as calories per mole 
of reactant, but from reference 39 it appears that the 
enthalpy change is for Eq. 1 as it is shown. 
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ture in question. Equation 3 gives the amount 
of heat, in calories, necessary to heat a mole 
of steam from a temperature 7, to T. 


AH, = 7.219 [T — T,] + 1.187 x 10° [T? - 79] 


+8.9x 10° [7° — T3]* cal/mole (3) 
Since Eq. 2 was derived from free-energy data, 
it is necessary to use the same standard state 
and 7; value in Eq. 3 as were usedin calculating 
Eq. 2. The authors use a value of 298°K for 
T;. The difference between the heat liberated, 
AHp, and the heat passing to the steam, AHz, 
is called the “heat entering the solid phase,”’ 
HESP. Since AWis the oxygen weight gain in 
time ¢ (grams of oxygen per square centimeter), 
the product of AW and HESP is the heat liber- 
ated per unit area if HESP is expressed in 
calories per gram of oxygen. 

The heat balance can be expressed as fol- 
lows: 


Specific heat = decay heat + oxidation heat 


or 
AT x Choy aK W =H) x tx W + HESP 
x [R(T) (t + to)” — R(T) 127] (4) 


In this equation: 

1. AT is the temperature rise of the element 
during time ¢, °K. 

2. Cphz, is the specific heat of the cladding, 
cal/g/°K. 

3. W is the weight of a half-section, g/cm’. 

4. Hp is the decay heat liberated, cal/g/sec. 

5. R(T) represents the total number of 
grams of oxygen reacted prior to the tempera- 
ture interval in question. 

6. ¢) is the time (in seconds) at the new tem- 
perature that it would take to oxidize the sec- 
tion to the extent it has reached prior to the 
temperature interval in question. 

7. R(T)(t + t,)"!7) is the total amount of oxi- 
dation that has taken place, including the tem- 
perature interval in question. 





* The authors of reference 38 noted that the units 
on this equation are not correct since an error in the 
original assumption resulted in the lower limit of the 
equation being taken for 373°K rather than 298°K. Ac- 
cordingly, the Lustman values calculated are slightly 
high at the lower temperatures. 


8. The difference between 7 and 5 is equal 
to the amount of oxidation taking place during 
the temperature interval in question. 


Equation 4 is solved for ¢ by trial and error, 
and the amount of oxidation taking place during 
the interval is calculated. This procedure is 
repeated until the melting temperature of the 
cladding is reached. At this point the left side 
of Eq. 4 is replaced by the heat of fusion, and 
the calculations are continued. The calcula- 
tions would continue until the core had melted 
or some external mechanism (e.g., water 
sprays) had terminated the excursion. 

The second part of the problem is started 
when the cladding ruptures exposing the fuel 
to the water. In the particular reactor design 
studied in reference 38, the fuel, an alloy of 
uranium, had a lower melting temperature than 
Zircaloy. Hence, when the Zircaloy melted, it 
was assumed that this molten fuel suddenly 
dropped into the pool of water. If the fuel were 
UO,, this part of the problem, of course, would 
not exist. In addition, molten drops of Zircaloy 
would be dropping off fuel elements and into 
the water pool at the bottom of the pressure 
vessel. The heat balance is again the basic 
equation to solve, using the method of dividing 
the total temperature interval into increments 
of width AT. It is assumed that heat is lost 
from the molten drops entirely by radiation, 
and the balance equation becomes 


Hrad =Hp + py + Hsy 


where H,.q = energy transferred by radiation 
Hp =energy from metal-water reac- 
tion 


Hpy = decay energy 
Hsy =energy liberated in cooling an 
amount AT 


Several additional complications arise, how- 
ever. The molten metals form drops, which 
must be sized, and the calculations require 
emissivities, which must be determined. Also, 
as the drop oxidizes, the surface area available 
for further oxidation decreases. Droplet sizes 
can be approximated from interfacial tension 
data, and the reduction of surface area can be 
calculated as the oxidation proceeds. 

When the calculation is complete, the amount 
of metal-water reaction is known, subject to 
the approximations used in the derivations and 
in data such as emissivities and interfacial 
tension. Tables VI-2 and VI-3, for example, 
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Table VI-2 TIME AND PERCENTAGE REACTION FOR ZIRCALOY-2-CLAD FUEL ELEMENT™® 
(Decay Heat = 0.374 cal/g/sec) 

















Total 
HESP, * oxidation, Total 
Temp., °K cal/g O k(T) n(T) to, sec t,sec g O/cm? oxidation, % 
800-1000 4090 $:3.% 16" 0.33 0 25.5 9.03 x 1074 0.45 
1000-1200 3920 4.15 x 10 0.38 42 29.3 1.64 x 1073 0.81 
1200—1400 3740 6.31 x 1074 0.462 7.9 20.1 2.95 x 107% 1.5 
1400—1600 3580 9.6 x 10 0.552 KY J 11.0 4.85 x 1073 2.4 
1600—1800 3410 1.52 x 107 0.66 6.1 4.6 7.26 x 107 3.6 
1800-2000. 3230 2.5 x 1073 0.784 3.9 2.0 10.0 x 107% 5.0 
2000—2200 3040 4.36 x 1079 0.945 2.4 0.7 12.7 x 1073 6.3 
2200— molten 3040 4.36 x 1073 0.945 3.1 3.0 24.04 x 107 11.9 
* HESP = heat entering solid phase; Lustman’s values. 
Table VI-3 COOLING DATA FOR A DROPLET OF ZIRCALOY-2™® 
(Radius = 0.215 cm, Emissivity = 0.67, Decay Heat = 0.374 cal/g/sec) 
HESP Total 
Temp., °K (cal/g O) k(T) n(T) t,sec Total g O reaction, % 
2200—2000* 3040 4.36 x 1073 0.945 4.0 9.38 x 1078 9.8 
2000—1800 3230 2.5 x 1073 0.784 1.1 10.84 x 107% 11.3 
1800-1600 3419 1.52 x 1073 0.660 1.4 11.85 x 107% 12.4 
1600— 1400 3580 9.6 x 10 0.552 2.1 12.62 x 1073 13.2 
1400—1200 3740 6.31 x 10~¢ 0.462 3.3 13.20 x 107° 13.8 
1200—1000 3920 4.15 x 1074 0.38 6.0 13.63 x 1079 14.2 
1000-800 4090 3.1 x 1074 0.33 13.8 14.6 





*Includes solidification. 


show the reaction history of a Zircaloy-2-clad 
element and a drop of molten Zircaloy-2. 


SRE Shutdown 


During a normal power run on the SRE in 
May, tetralin coolant from the main primary 
pump shaft seal leaked into the primary sodium 
system. The reactor was shut down and the 
pump seal coolant system converted from tetra- 
lin to NaK. It was started up again in mid-July 
to remove tetralin decomposition products 
through careful low-power operation and cold 
trapping. Anomalies were noted during the run 
as expected. When several weeks of operation 
failed to clear the anomalies, the reactor was 
shut down to inspect fuel elements. 

Since shutdown in late July, all the fuel ele- 
ments have been removed from the reactor. 
Twelve of the 43 were found to be severely 
damaged by high temperature. The mechanism 
of damage has been established as partial 
blockage of reactor coolant flow passages by 
tetralin decomposition products, followed by 
local overheating of the fuel elements, resulting 
in formation of an iron-uranium eutectic (m.p. 
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14.02 x 1078 


1340°F) and parting of the type 304 stainless- 
steel cladding. 

Samples of the fuel, cladding, and other com- 
ponents of the fuel assemblies have been metal- 
lurgically examined, and coolant and coolant 
contaminants have been chemically analyzed. 

The results of these investigations will be 
given in a report by the SRE Ad Hoc Review 
Committee scheduled for publication in Decem- 
ber. The report will be reviewed in the March 
1960 issues of Power Reactor Technology and 
Nuclear Safety. 
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SHIELDING 





Review of Shielding Art 


In reference 1 is a review of analytical shield- 
ing methods. Each method is discussed, and 
some are used in predicting the results of 
experiments. The comparisons cover multi- 
region, multilayer, and multimaterial shields; 
ducts and porous plugs; air and ground scatter- 
ing; shield heating; and other subjects of lesser 
interest in power reactor shield design. All 
comparisons include neutron and gamma dose 
rates with the addition of two comparisons of 
measured and calculated fast-neutron energy 
spectra. The calculational methods are not 
described in detail, but all methods are ref- 
erenced. 
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Although the specific configurations investi- 
gated are of little interest to the designer of 
large power reactors, the calculation techniques 
used give striking agreement with experimental 
results, even in situations which are quite com- 
plex. There is little doubt that these techniques 
would be valuable in the civilian power reactor 
program and would constitute a real advance 
over those presently in use. 
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Radiation Effects . 
on Semiconductors 


Nuclear radiation poses a severe problem inthe 
application of semiconductors to control andin- 
strumentation systems in the nuclear field. Be- 
cause of the susceptibility of semiconductor 
materials to radiation damage, an understanding 
of the permissible limits of radiation exposure 
and of the changes in characteristics of the 
semiconductor is required. The Radiation Ef- 
fects Information Center (REIC) has reportedin 
two earlier memorandums’ on some of the ef- 
fects of radiation on semiconductor materials 
and more recently has issued a memorandum 
on the effects of radiation on semiconductor 
diodes.’ Transient, semipermanent, and perma- 
nent diode parameter changes occur as a result 
of exposure. 

Reference 3 has made the following summary 
of semiconductor diode behavior under nuclear 
radiation: 


1. Fast neutrons affect semiconductor diode for- 
ward characteristics through changes in the bulk 
material, and they affect the reverse characteris- 
tics through changes in the junction region. 

2. Where fast- and slow-neutron fluxes are of 
the same order of magnitude, the slow neutrons are 
not a major factor contributing to the deterioration 
of semiconductor devices. 

3. Gamma radiation causes transient surface 
damage in semiconductor diodes which results in 
changes in leakage currents; recovery may occur 
when the device is removed from the radiation field. 

4, Germanium diodes fail in the reverse direc- 
tion because of changes in reverse resistance of 
several orders of magnitude after exposure to inte- 
grated neutron fluxes of 10'3 fast neutrons/cm? or 
higher. 

5. If germanium diodes are employed in elec- 
tronic circuitry exposed to integrated fluxes as 
great as 10'* fast neutrons/cm?, the circuit design 
must be such that (a) thecircu.t performance is not 
seriously affected by changes in the reverse char- 
acteristics of the germanium diodes and (b) the 
diodes are protected individually from self- 
destruction resulting from excessive power dissipa- 


tion when the radiation deteriorates the reverse 
characteristics. 

6. Failure in silicon rectifying devices results 
primarily from deterioration of the forward char- 
acteristics. For applications where an increase in 
forward resistance does not seriously affect opera- 
tion, silicondiodes probably can be used to 10" fast 
neutrons/cm?, Equipment which utilizes silicon di- 
odes and which will be exposed to integrated fluxes 
to 10'3 fast neutrons/cm? must be capable of per- 
forming satisfactorily with diode reverse leakage 
currents as large as approximately 0.2 wa at about 
30 volts. 

Some silicon diode types are highly radiation re- 
sistant and have usable forward characteristics to 
10'3 fast neutrons/cm? in addition tousable reverse 
characteristics. If the forward characteristics are 
important to equipment performance, radiation- 
resistant units such as the NSP type 1N137A silicon 
diode should be employed. 

7. Radiation-invoked ionization of air molecules 
adjacent to semiconductor diodes will result in col- 
lection of ionized molecules on bare diode leads, 
other conduction surfaces, and terminals and will 
cause currents to 2 ordersof magnitude larger than 
those created inside the diode envelope by radiation. 
Proper insulation virtually eliminates such cur- 
rents. 

8. A major factor determining the radiation re- 
sistance of semiconductor diodes appears to be the 
manufacturing process employed. 

9. Materials currently used in the fabrication 
of diodes are relatively. ineffectual for shielding 
against the effects of neutron and gamma radiation. 


Radiation Effects 
on Magnetic Materials 


Magnetic materials are used in motors, am- 
plifiers, transducers, and other actuators and 
instruments in the nuclear field. When these 
materials are exposed to large amounts of neu- 
tron radiation, some changes in properties of 
the materials can be expected. 

The REIC has summarized in a technical 
memorandum‘! some of the effects of nuclear 
radiation on the structure-sensitive properties 
of magnetic materials, such as permeability, 
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remanence, and coercive force. Radiation ef- 
fects on magnetic core materials of nickel-iron, 
silicon-iron, aluminum-iron, and ferrite were 
reviewed. Important changes in the magnetic 
properties of the nickel-iron alloys were noted 
after exposure to 10'” fast neutrons/cm’. No 
damage was observed in the operation of mag- 
netic amplifiers to exposures of 10'° fast 
neutrons/cm’. 


Corrosion of Stainless Steel 
in a Nitrogen Atmosphere 


Battelle Memorial Institute (BMI) has recently 
reported on corrosion tests on stainless-steel 
fuel-element cladding materials for the Gas- 
cooled Reactor Experiment (GCRE).° The fuel 
elements for this reactor experiment will be 
cooled by nitrogen, and the clad surface tem- 
perature is expected to be as high as 1650°F. 
Preliminary screening of stainless steels having 
satisfactory high-temperature properties were 
conducted at atmospheric pressure andelevated 
temperature with nitrogen gas containing addi- 
tions of oxygen and water vapor. The results of 
these tests are shownin Tables VIII-1 and VIII-2. 

As a result of the screening tests, AISI type 
318 stainless steel was selected for further 
study. A general indication of the corrosion 
rate of this material as a function of tempera- 
ture and time is provided by the data of Tables 
VIII-3 and VIII-4. The addition of small amounts 
of oxygen was shown to inhibit nitriding and result 
in the preferential oxidation of the stainless 
steel. 

The data obtained are encouraging for short- 
time use of this material at high temperatures. 


Table VIII-1 


However, at 1500°F, the preliminary data indi- 
cate that this material is still not suitable for 
long-term service in a power reactor using 
nitrogen coolant unless the corrosion rate drops 
significantly with time. 

The general conclusions made by BMI as a 
result of the screening study are quoted below.° 


1. Type 318 stainless steel has better resistance 
to oxidizing attack in wet and dry oxygen-containing 
nitrogen atmospheres than do type 347, type 316, 
Timken 16-15-6, and Timken 16-25-6 alloys. Type 
347 stainless was nearly as satisfactory as type 
318. 

2. The relationship between depth of attack and 
time of exposure seems to be diffusion controlled, 
indicating that the rate of attack decreases with 
time. This relationship may be disturbed adversely 
by spalling, the occurrence of which would cause 
the rate to increase. 

3. Within wide limits the surface finish does not 
affect the rate of oxidation of type 318 stainless in 
wet ordry oxygen-containing nitrogen atmospheres. 

4. Nitriding reaction is inhibited when small 
amounts of oxygen are present in a nitrogen atmos- 
phere. Some nitrogen diffuses into the metal, but 
nitrides do not appear to form. The resulting oxi- 
dation attack is much less severe than the nitriding 
attack. 

5. The concentration of oxygen in nitrogen which 
was found to be optimum from the standpoint of 
minimum oxidation attack and inhibition of nitriding 
was 0.5 vol.%. The concentration of oxygen found to 
be most effective was 5 vol.%, followed by 0.1 and 
0.05 vol.%. The results obtained when 5 vol.% oxy- 
gen was added to the atmosphere would warrant 
trying out dry and moist air atmospheres in the 
future. 

6. The addition of 0.6 vol.% water vapor to a dry 
nitrogen atmosphere does not affect to any great 
extent the rate of attack of type 318 stainless as 


MAXIMUM DEPTH OF OXIDATION OF VARIOUS FERROUS ALLOYS EXPOSED FOR 


100 HR IN NITROGEN CONTAINING SMALL AMOUNTS OF OXIDIZING GASES*° 





Maximum depth of penetration after 100 hr in nitrogen having addition shown (mils) 











0.2 vol.% 0.05 vol.% 
0.07 vol.% H,O H,O oxygen 0.1 vol.% oxygen 0.5 vol.% oxygen 
Alloy 1500°F 1650°F 1800°F 1650°F 1650°F 1800°F 1500°F 1650°F 1800°F 1500°F 1650°F 1800°F 

Type 316 S.S. t 0.5 1.6 t 0.2 1.0 0.05 t 2.0 0.1 0.15 t 
Type 318 S.S. 0.05 0.05 5.0 0.15 0.05 0.2 0.05 0.05 0.25 0.05 0.05 0.15 
Type 347 S.S. Tt 0.05 5.0 sf 0.6 0.4 0.05 Tt 0.3 0.05 0.2 t 
Timken 

16-15-6 0.3 0.1 1.3 0.5 0.15 0.9 0.1 0.3 1.0 0.2 0.6 0.9 
Timken 

16-25-6 0.3 0.3 0.6 0.3 0.3 0.5 0.15 0.3 0.5 0.1 0.4 0.6 





*Specimen surface was ground to 240-grit finish. 
+ Material was not tested. 
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Table VIII-2 MAXIMUM DEPTH OF OXIDATION OF 
VARIOUS FERROUS ALLOYS EXPOSED TO OXIDIZING 
NITROGEN ATMOSPHERES FOR 500 HR AT 1650°F** 





Maximum depth of penetration 
in alloy shown (mils) 








Atmosphere 
composition, Type Type Type Timken Timken 
vol.% 316 318 347 16-15-6 16-25-6 
Nitrogen 
+0.07 H,O 1.3 0.6 1.5 0.6 0.7 
Nitrogen 
+0.1 oxygen 1.2 0.6 1.0 0.6 0.7 





*Specimens were ground to 240-grit finish and annealed 
at 2150°F in hydrogen for 1% hr. 


Table VOI-3 MAXIMUM DEPTH OF ATTACK OF 
TYPE 318 STAINLESS STEEL EXPOSED AT SEVERAL 
TEMPERATURES TO DRY NITROGEN ATMOSPHERES 

CONTAINING 0.5 VOL.% OXYGEN® 








Preanneal Exposure Depth of 
Speci- temp.,* temp., Exposure attack, 
Test men ri 2 ss time, hr mils 
1 5 2150 1750 215 0.45 
6 2150 1750 376 0.7 
44 2150 1750 479 0.7 
7 2150 1750 580 0.95 
41 2150 1750 683 1.1 
27 2150 1750 843 1.3 
8 2150 1750 1059 21 
2 17 2150 1650 216 0.4 
18 2150 1650 478 0.5 
37 2150 1650 980 0.75 
76 2150 1650 1032 1.6 
58 2150 1650 2468 1.0 
20 2150 1650 2626 0.8 
40 2150 1650 3180 1.8 
19 2150 1650 3658 2.9 
3 54 2450 1650 238 0.7 
77 2450 1650 417 0.5 
55 2450 1650 662 0.9 
73 2450 1650 1089 0.8 
56 2450 1650 1334 1.0 
4 64 2450 1500 258 Nil 
80 2450 1500 417 0.25 
65 2450 1500 930 0.2 
72 2450 1500 1176 0.6 
66 2450 1500 1347 0.4 
67 2450 1500 1434 0.5 





*All specimens were annealed in hydrogen. Those an- 
nealed at 2150°F were held for 2 hr and furnace-cooled. 
Specimens annealed at 2450°F were held for 5 hr, cooled to 
1850°F, and held for 2 hr, and then furnace-cooled. 


measured by the maximum depth of penetration, but 
it does result in a more nonuniform attack. 

7. On the basis of limited work, it can be said 
that nitrogen atmospheres with addition of small 
amounts of carbon monoxide cause shallower oxida- 


Table VIII-4 MAXIMUM DEPTH OF ATTACK OF TYPE 318 
STAINLESS STEEL EXPOSED AT SEVERAL TEMPERA- 
TURES TO NITROGEN ATMOSPHERES CONTAINING 
0.5 VOL.% OXYGEN PLUS 0.6 VOL.4 WATER VAPOR? 











Preanneal Exposure Depth of 
Speci- temp.,* temp., Exposure ' attack, 
Test men = °F time, hr mils 
5 29 2150 1800 50 0.6 
30 2150 1800 210 2.8 
31 2150 1800 616 5.4 
32 2150 1800 808 5.6 
6 33 2150 1650 50 0.15 
34 2150 1650 210 0.4 
79 2150 1650 417 0.9 
35 2150 1650 616 1.1 
74 2150 1650 1080 0.7 
36 2150 1650 1498 0.75 
52 2150 1650 1971 1.8 
50 2150 1650 2120 2.0 
43 2150 1650 2377 1.8 
7 84 2450 1500 65 Nil 
68 2450 1500 210 Nil 
83 2450 1500 319 9.3 
*All specimens were annealed in hydrogen. Those an- 


nealed at 2150°F were held for 2 hr and furnace-cooled. 
Specimens annealed at 2450°F were held for 5 hr, cooled to 
1850°F, held for 2 hr, and then furnace-cooled. 


tion attack than do nitrogen atmospheres containing 
small amounts of oxygen. 

8. Several conclusions concerning the mechanism 
of attack can be drawn from the X-ray diffraction 
and fluorescence analyses of the oxidation products. 
These are: 

a. The initial oxidation of the alloy produces an 
oxide film in which the major phase is of the R,O; 
type. An- oxide of the R,O, spinel type is a minor 
component in the initial film. The metal constitu- 
ents of the initial film are presumed to be essen- 
tially those of the alloy. 

b. Preferential diffusion of chromium and man- 
ganese from the base metal into the oxide layer 
during the next step results in the formation of in- 
creasing amounts of the R,O, type oxide. 

c. The amount of preferential diffusion is reduced 
significantly by the addition of water vapor to the 
oxidizing gas mixtures. Diffusion is only slightly 
decreased by increasing the oxygen content of the 
gas over the range 0.05 to 5.0 vol.%. 

d. At 1650°F, the composition of the oxide layer 
reaches a maximum in chromium and manganese 
and a minimum in iron and nickel after 500 to 600 


hr of oxidation. 

e. The composition of the adherent oxide, after 
500 to 600 hr of oxidation, varies considerably due 
to spalling of parts of the oxide layer. 

f. Oxidation of the thicker noncoherent oxide 
layer results in the formation of a chromium-rich 
R,O, oxide. 
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9. Specimens that were preannealed at 2150°F in 
hydrogen appeared to be slightly more resistant to 
attack than specimens preannealed in hydrogen at 
2450°F. 
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LIGHT-WATER REACTORS 





The Belgian BR-3 


The Belgian BR-3 nuclear power plant is pres- 
ently under construction at Mol, Belgium, and 
should commence initial operation during the 
summer of 1960. Reference 1 contains recent 
information on the design of the plant. The 
following organizations are responsible for the 
designated phases of the operation. 


Centre d’Etude de 1’Energie Nu- 
cléaire (CEN), Brussels, Bel- 


Plant owner 


gium 
CEN’s consulting Bureau d’Etudes Nucléaires 
engineers (BEN) 


Société Belge pour 1’Industrie 
Nucléaire (Belgo-Nucléaire) 
Société de Traction et d’Elec- 

tricité (TE), Brussels, Belgium 
Westinghouse Electric Corp., 
clear equipment Atomic Power Division 
and systems (WAPD), Forest Hills, Pa. 
Architect-Engineer Gibbs and Hill, New York, N. Y. 
for Westinghouse 
Plant container 
designer 


Designer of nu- 


Sargent & Lundy, Chicago, Ill. 


The reactor is a PWR type employing light 
water as the coolant and moderator. Other 
design details can be found in Table IX-1. 
Boric acid solutions are used for two purposes 
in the design. The first system injects boric 
acid into the primary coolant to achieve shut- 
down to room temperature for refueling pur- 
poses, since the control rods are not adequate 
to keep the clean, cold reactor in a subcritical 
condition. A second system supplies borated 
water to the reactor for core cooling in the 
event of loss of coolant through a rupture in 


Table IX-1 


DESIGN DETAILS FOR BR-3 REACTOR 





Reactor thermal power 

Gross electric power 

Net electric power 

Steam pressure at throttle 
Steam temperature at throttle 
Steam flow rate 

Feed-water temperature 


Reactor coolant conditions (nominal): 


Pressure 
Inlet temperature 
Outlet temperature 
Flow rate 
Average heat flux 
Peak heat flux 
Fuel-element geometry 


No. of tubes per assembly 
No. of assemblies per core 
Enrichment: 

Outer zone 

Inner zone 
Core life 


Control 


Reactivity change; cold, clean to 
hot, end of life 

Containment-shell dimensions, H/D 

Containment-shell maximum design 
pressure 


40,900 kw 
11,500 kw 
10,500 kw 
525 psia 
472°F 
154,500 lb/hr 
330°F 


2000 psia 

491°F 

515°F 

5 x 10° lb/hr 

94,200 Btu/(hr)(sq ft) 

441,800 Btu/(hr)(sq ft) 

Stainless-steel tubes 
containing UO, 
pellets 

~110 

32 


4.4 wt.% 

3.7 wt.% 

7000 equivalent full- 
power hours 

12 rods plus boric 
acid 

19 per cent 


107 by 54 ft 
45 psig 





the primary-coolant plumbing. This latter sys- 


tem 


is actuated manually when the system 


pressure has fallen to 800 psia. 


Reference 


1. P. Dozinel and J. Storrer, Safety Aspects of the 
Belgian BR-3 Reactor, presented at the Interna- 
tional Symposium on the Safety and Location of 


Nuclear Plants, June 1959. 
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One of the most important aspects of boiling- 
water reactors is the stability of power output. 
Considerable progress has been made in this 
field; the subject is reasonably well understood 
qualitatively, and it is possible to design large 
high-output boiling reactors, as well as the 
smaller reactors, without fear of encountering 
instability in the operating range. Nevertheless, 
the designer is still chary of those design fea- 
tures which are thought to promote instability, 
and, as a result, the design areais perhaps more 
restricted than is necessary. 


The subject of boiling-reactor stability has 
been discussed in previous issues of this Review 
(Power Reactor Technology, 1(1): 14-16 (Decem- 
ber 1957); 1(3): 21-22 (June 1958); 1(4): 10-14 
(September 1958)). Recently, several new re- 
ports have been issued which make contributions 
to the understanding of the subject. 


Design Considerations 


The closest relation between questions of sta- 
bility and reactor design may be expressed in 
terms of the steam coefficient of reactivity, 
sometimes called the void coefficient of reac- 
tivity. When steam exists in a reactor, the 
reactivity will, in general, differ from the 
reactivity when there is no steam present. If the 
steam in the reactor core occupies a fraction 
AV/V of the volume available for occupancy 
by H,O and, in doing so, causes a change of 
reactivity Ak or:/k es, then the steam coefficient 
of reactivity is defined as (Ak ef¢/keft)/(AV/V). 
If the coefficient is negative, as it must certainly 
be if the reactor is to be inherently stable, then 
any increase in reactor power which will in- 
crease the steam content of the reactor will 
decrease the reactivity; this, in turn, will tend 
to decrease power. In other words, the reactor 
will tend to regulate its power to a constant 
value if no reactivity changes are imposed by 
outside sources. The regulation will not occur 
instantaneously since there are various thermal 
and hydraulic time delays in the process. 
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It is both customary and useful to consider 
(1) the externally imposed reactivity variations 
to constitute the input signal, and (2) the reactor 
power output to constitute the response, of a 
servo loop which includes those reactor charac- 
teristics involved in the reactivity-power rela- 
tion (Fig. 12). In this representation the net 
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Figure 12—Representation of the reactivity-power 
relation in terms of a servo system. 


reactivity is applied as a driving force to the 
neutron chain reaction, producing the appro- 
priate variation in neutron flux and, hence, in 
reactor power. Those processes (e.g., change of 
steam content) which, as a result of a power 
change, produce a reactivity change are con- 
sidered to constitute a feedback circuit, whose 
output, the feedback reactivity, combines alge- 
braically with the externally applied reactivity 
(the input signal); this combination yields the 
net reactivity applied to the neutron chain 
reaction. In this representation the thermal and 
hydraulic time delays appear as phase lags in 
the frequency characteristics of the feedback 
loop. It is characteristic of such systems that 
they may oscillate unstably for certain combina- 
tions of phase shifts and gain of the feedback 
circuit. In general, the tendency toward in- 
stability will increase as the gain is increased. 
In the case of any particular boiling reactor, 
the steam content of the reactor is roughly pro- 
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portional to the operating power. Consequently, 
the fractional change in reactivity which results 
from a fractional change of power (which may 
be considered the gain of the feedback loop) is 
roughly proportional to the product of the power 
times the steam coefficient of reactivity. That 
is to say, if time lags in the system are such 
that instability is possible, then the point of 
instability will be approached more and more 
closely as reactor power is increased. Since, in 
general, the designer wishes to achieve a rela- 
tively high power density in the reactor, and 
since he is aware that a stability limit probably 
exists which acts as a limit on the product 
(power) x (steam coefficient of reactivity), his 
incentive is to make the steam coefficient of 
reactivity as low as possible. Hence there is a 
tendency to design boiling reactors for a steam 
coefficient of reactivity whichis low in magnitude 
and negative in sign. 

If one considers the reactivity change which 
may occur in a water-moderated reactor when 
some of the water is displaced by steam —the 
reactivity change which determines the steam 
coefficient of reactivity — several contributions 
to the change may be recognized. The parasitic 
neutron absorption by water will be decreased, 
and the fast-fission factor may be slightly in- 
creased; both of these effects make a positive 
contribution to the reactivity change. The reso- 
nance absorption in U?** will increase, and the 
leakage of neutrons out of the reactor and the 
leakage of neutrons into control rods and other 
lumped absorbers which may be present will 
increase; these effects make negative contribu- 
tions to the reactivity change. The relative 
importance of these effects will vary with the 
characteristics of the reactor design. 

In particular, if the reactor fuel contains a 
fertile material, the relative importance of 
parasitic absorption by water and of resonance 
absorption by the fertile material can be ad- 
justed by adjusting the ratio of fuel to water. 
Higher fuel-to-water ratios will accentuate the 
effect of resonance absorption and cause the 
negative component of the steam coefficient to 
predominate, whereas low fuel-to-water ratios 
will accentuate the parasitic absorption of water 
and cause the positive component to predominate. 
In such a reactor the net steam coefficient of 
reactivity can be adjusted to any desired value 
over a rather wide range by suitable adjustment 
of the fuel-to-water ratio. It is evident from the 
same consideration that the steam coefficient of 


reactivity will not remain constant as the amount 
of steam in the reactor increases, but that it 
will become more and more strongly negative as 
the average density of water in the reactor is 
decreased by an increase in steam content. A 
similar effect will result from a temperature 
increase which decreases the water density. In 
a high-pressure boiling reactor the decrease in 
average water density may be some 25 per cent 
as a result of the temperature change from room 
temperature to operating temperature, and a 
further decrease of something like 15 per ceni 
may result from the operating steam content. 
Consequently, the steam coefficient may be much 
more strongly negative under full-power operat- 
ing conditions than it is at room temperature. 
Characteristically, the reactor designer is faced 
with a choice between (1) designing for a very 
low negative coefficient at room temperature anc 
accepting a somewhat higher negative coefficient 
than he desires under operating conditions and 
(2) designing for the desired lownegative coeffi- 
cient under operating conditions and accepting a 
positive coefficient at room temperature. If the 
latter course is taken, a hazard question may 
arise because of the possibility of an autocata- 
lytic reactivity accident at lowtemperature, and 
it may be necessary to forego operation of the 
reactor at low temperature and provide auxiliary 
heating to bring the reactor up to an acceptable 
operating temperature. 

In so far as the fuel-to-water ratio must be 
chosen to yield the desired steam coefficient of 
reactivity, other design goals which are affected 
by fuel-to-water ratio may be compromised. 
One of the important goals of this kind is high 
conversion ratio. If there is no restriction on 
fuel-to-water ratio, the conversion ratio may 
be increased by increasing the fuel-to-water 
ratio. The improvement in conversion ratio re- 
sults from increased resonance absorption in the 
fertile material, which, in turn, allows the use 
of a higher fuel enrichment with correspondingly 
smaller fractional loss of neutrons to parasitic 
absorbers. Fortunately, the fuel-to-water ratios 
which yield acceptable steam coefficients of 
reactivity also give reasonably high conversion 
ratios (in the neighborhood of 0.7) if low- 
absorption (Zircaloy) jackets are used for the 
fuel. Nevertheless, greater freedom in adjust- 
ment of conversion ratio through adjustment of 
fuel-to-water ratio (or more precisely through 
adjustment of the resonance escape probability) 
would be desirable. 











44 PROGRESS ON SPECIFIC REACTOR TYPES 


It must be evident, from the above considera- 
tions, that overconservatism in the physics de- 
sign of the boiling reactor, i.e., the striving for 
an extremely low steam coefficient of reactivity 
in order to avoid encroachment on the gray area 
in which stability is questionable, entails penal- 
ties in reactor cost and performance. Equally 
important penalties may arise from overcon- 
servatism in the thermal design. 

The conservatism in the design of the Experi- 
mental Boiling Water Reactor (EBWR) has been 
thoroughly demonstrated, and such a degree of 
conservatism was certainly appropriate for an 
early boiling reactor built as an experiment. It 
is not appropriate for a commercial power 
reactor. Such a reactor is ordinarily built to 
a specified power capacity, and for economic 
reasons it should be the least reactor which will 
yield that capacity reliably and safely. Further, 
in such a practical machine, other factors which 
limit power output will, in the economic design 
of the reactor, be harmonized with the supposed 
stability limitations and will, thereafter, prob- 
ably limit the reactor capacity to the design 
value. For example, in order to minimize fuel 
cost, the UO, fuel pellets will probably be made 
as large as is permissible within the limitations 
of center melting at the design power density. 
Thereafter the fuel temperature will limit the 
power output even if operation of the reactor 
demonstrates that a considerably higher output 
would be permissible from the standpoint of 
stability. 

In summary, although boiling reactors can be 
designed on the basis of present knowledge to 
compare favorably with other reactor types, 
the lack of complete understanding of stability 
does prevent design optimization. The deficiency 
affects not only the design of specific reactors 
but also the forecasting of the ultimate value of 
the reactor type; it also affects the choices 
between types of boiling reactors (e.g., forced 
circulation versus natural circulation) for a 
given application. For these reasons, new in- 
formation on boiling-reactor stability is always 
of considerable interest to the reactor designer. 


Experimental Boiling Water 
Reactor 
References 1 to 3 are a series of reports 


related to the EBWR and tothe transfer-function 
measurements which were made on that reactor. 


Reference 1 describes the measurements them- 
selves; reference 2 describes the major dynamic 
characteristics of an EBWR type reactor from 
the theoretical point of view and derives ap- 
proximate expressions for the transfer functions 
of the important contributors to the dynamic 
behavior. In reference 3 the theoretical approach 
of reference 2, with some modification, is ap- 
plied to an analysis of the measured EBWR 
transfer functions; by asemiempirical approach, 
rational values for the major dynamic com- 
ponents are derived. 

The transfer function of a reactor is simply 
the frequency response of the reactor power, 
with respect to amplitude and frequency, to a 
sinusoidal oscillation of reactivity, whichis pro- 
duced in this case by the oscillation of a control 
rod. As is well known from other applications of 
frequency-response techniques, such measure- 
ments can be used to indicate the frequency 
regions of potential instability of a dynamic 
system when the system is operating at a gain 
lower than that which marks the threshold of 
instability. Such analyses can, of course, also 
be used as a basis for inferring the more 
fundamental dynamic characteristics of the 
system. 

As has been implied earlier, the dynamic 
response of a reactor to an impressed reac- 
tivity oscillation can be considered to be made 
up of two major components. The first of these 
is simply the response of the neutron chair 
reaction to the oscillation of reactivity. For 
small oscillations to appear in the frequency 
range of practical interest, it depends only on 
the delayed-neutron characteristics and the 
effective prompt-neutron lifetime. All other 
dynamic effects result from reactivity changes, 
which themselves are caused by the changes in 
reactor power, and add to or subtract from the 
impressed reactivity oscillation. These effects 
are properly thought of as feedback effects that 
may modify the reactor power response drasti- 
cally from the response of the “zero-power” 
neutron chain reaction. 

The relation between the chain-reaction re- 
sponse and the feedback effects can be expressed 
diagrammatically by the usual conventions of 
servomechanism theory as in Fig. 12. In this 
diagram, G represents the complex response, 
or the transfer function, of the neutron chain 
reaction, and H represents the net power- 
reactivity transfer function. That is to say, His 
the transfer function which shows the combined 
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effect of all the processes that may result in a 
change of reactivity because of a change in re- 
actor power. In general, the gain of all these 
feedback effects is proportional to the operating 
power level. Consequently, since the kinetic 
behavior of the neutron chain reaction is in- 
dependent of power level, it is a relatively 
simple matter to determine the G function 
explicitly simply by measuring the reactor 
transfer function at such sufficiently low power 
that the feedback effects are negligible. For 
this reason the G function is often referred to 
as the “zero-power” transfer function of the 
reaction. If G has been so determined, and if a 
further transfer-function determination is made 
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Figure 13—-EBWR transfer functions, 20 Mw(t), 300 psig. 


at an appreciable power level, then it is theoreti- 
cally possible to determine the H function for 
that power level by separating the G function 
from the reactor transfer function. The equation 
which relates the reactor transfer function, P, 
to G and H is 


— 
P "i+ GH 
It is to be understood that this is a symbolic 
equation; since each point of eachtransfer func- 
tion (at a given frequency) is specified in terms 
of an amplitude and a phase relation, the indi- 
cated mathematical manipulations are to be 
carried out in the complex plane. 
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, feedback transfer function (H). 


----, zero-power transfer function (G). —-—, power transfer function (P). Data are taken from 


reference 3. 
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Figure 14—-Bode diagram for EBWR feedback func- 
tion, 20 Mw(t), 300 psig. A, reduction of apparent 
steam coefficient of reactivity by “void flashing” and 
change of nonboiling length. B, reduction of apparent 
steam coefficient of reactivity by “void flashing” only. 
C, w=1/T,, where T, is a fuel heat-transfer time 
constant. D, w=1/T»,, where Tp, is a fuel heat- 
transfer time constant. E, w=1/7T, where 7 is the 
steam transit time constant.’ 


Needless to say, the determination of H by the 
above procedure requires rather accurate ex- 
perimental data, for in some cases the value 
of H results from rather small differences of 
rather large numbers. Reference 3 discusses 
this problem as well as the method used for 
smoothing the experimental data. The results, 
although they may have some imperfections, 
are extremely useful since they relate directly 
to those dynamic factors which determine the 
reactivity feedback, and they allow one to gaina 
considerable insight into the details of the dy- 
namic behavior and, to some extent, to test 
theoretical predictions of the individual factors. 
The results of these procedures are illustrated 
by Figs. 13 and 14 (from reference 3). 

The curves on the right in Fig. 13 are the 
experimentally determined reactor power trans- 
fer function. The plotted points are those ex- 
perimentally determined. On the left in Fig. 13, 
the derived feedback function H is shown. 
Again, the experimental points are marked by 
circles. On these same grids are plotted, as 
dashed lines, the zero-power or chain-reaction 
transfer function G. The solid lines of the H 


function are derived from a semiempirical 
theory, in which specific physical components 
of the feedback dynamics are recognized, and 
the form of their behavior is calculated from 
first principles; but, the particular constants 
which enter the description of their behavior 
are chosen to yield a good fit with the experi- 
mental H function. The lines plotted for the 
reactor transfer function in Fig. 13 are those 
derived from the solid H curves. 

An understanding of the effects taken intc 
account in the derivation of the H function can 
be gained from Fig. 14. In Fig. 14, the amplitude 
of the H function has been replotted along witl 
a straight-line approximation to the amplitude 
curve generally referred to as a Bode diagram. 
The straight-line approximation has the virtue 
of identifying certain individual feedback com- 
ponents and of marking on the frequency plot the 
frequencies which correspond to the inverse 
time constants characterizing those components. 
Starting from the left side of the figure and 
progressing tuward higher frequencies, the 
components of the feedback which have been 
taken into account can be identified. First, note 
the dashed line marked K5,V. Ko5pis the steam 
coefficient of reactivity expressed in dollars of 
reactivity per cubic foot of steam in the reactor 
core. V is called the power void coefficient. It 
is expressed as the steam content of the reactor 
and is expressed in cubic feet per megawatt of 
reactor power, multiplied by the operating power 
level in megawatts to give a factor consistent 
with the power-independent analysis. The prod- 
uct K5,V is simply the steady-state change in 
reactivity which would be produced by unit 
fractional change in the power level. Ks, maybe 
considered the steady-state steam coefficient of 
reactivity. Under the dynamic condition, how- 
ever, the effect of this coefficient is very much 
reduced, as indicated by the difference between 
the dashed line and the left-hand portion of the 
straight-line approximation. This reduction re- 
suits from two effects, both of which have to do 
with the change in reactor pressure with power 
level. As the pressure increases because of a 
power increase, the saturation temperature in- 
creases also; in addition, some of the steam 
already formed in the fuel channels must con- 
dense to raise the temperature of the water in 
the channels to the higher boiling point. This 
effect is referred to in reference 4 and in 
Fig. 14 as the “void flashing” effect. The sec- 
ond effect is a change in the nonboiling length 
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in the coolant channel. As the pressure in- 
creases, the subcooling of the water entering 
the coolant channels will be increased, and the 
average length of the nonboiling region in the 
fuel channels will increase appropriately. Con- 
versely, if power were reduced, the pressure 
would fall, there would be some flashing of water 
to steam in the fuel channels, and the average 
length of the nonboiling region in the fuel channels 
would decrease. Both of these effects are 
counter to that of the primary steam coefficient 
of reactivity. If the power change in question 
were a step change which remained in effect 
indefinitely, then the pressure would, in the 
case postulated, eventually reach an equilibrium 
value, the total body of water would again reach 
thermodynamic equilibrium with the steam, and 
the effects on the steam content of the reactor 
would disappear except for a small residual 
effect on the length of the nonboiling region 
caused by the slight change in subcooling of the 
reactor feed water. That is to say, the effects 
discussed are not directly related to pressure 
per se, but they are relatedtothe rate of change 
of pressure. 

The effect on nonboiling length depends upon 
the introduction into the fuel channels of the 
subcooled recirculating water. If the frequency 
of oscillation is increased to a point at which 
the sinusoidal period is comparable with the 
recirculation time of the water, then the effect 
begins to decrease. This decrease is represented 
on the straight-line approximation by the line 
segment covering the frequency range from about 
1.3 to 3.0 radians/sec. At higher frequencies the 
effect of variation of the nonboiling length has 
disappeared in the straight-line approximation, 
and the feedback magnitude is almost as great 
as Ks5p, x V, being smaller only because of the 
residual steam flashing effect. This effect does 
not depend on the water recirculation but de- 
pends only upon the rate of change of steam 
pressure in addition to those variables which 
affect steam generation. The effect therefore 
persists to higher frequencies. 

At frequencies of about 4.5 and7 radians/sec, 
two time constants, 7; and Ty,, enter which 
introduce phase lags in the formation of steam. 
These two time constants include the effects of 


heat storage in the fuel plates as determined by 
the heat capacity in the plates and the heat flow 
characteristics which are determined both by the 
thermal conductivity of the plates and cladding 
and by the heat-transfer coefficient from fuel 


plate to water. However, the effects of thermal 
conductivity and heat-transfer coefficient are 
not identified explicitly. The thermal inertia 
identified by these time constants reduces the 
magnitude of the feedback function rapidly, as 
the frequency is increased. At a frequency of 
about 17 radians/sec, a further time constant, 
T, enters as a further phase lag. 7 is identified 
as the time constant related to the steam transit 
time in the reactor, that is, the time associated 
with the escape of generated steam from the 
reactor. At frequencies higher than about 20 
radians/sec, the combined effects of thermal 
time lags and steam transit time have reduced 
the amplitude of the feedback function to a 
practically negligible point; at higher fre- 
quencies, the reactor transfer function (the 
right-hand side of Fig. 13) would be expected 
to be identical with the zero-power transfer 
function to within the limits of observation. For 
frequencies greater than about 20 radians/sec, 
the self-regulating characteristic of this boiling 
reactor is virtually absent. 

Analyses similar to that represented by Fig. 13 
are reported in reference 3 for a number of 
reactor operating conditions, and the time con- 
stants and associated feedback parameters are 
derived for the several operating conditions by 
making experimental fits to the feedback transfer 
functions. The data so derived are of consider- 
able interest and should provide further insight 
into the behavior of the boiling reactor, even 
though a rather simplified picture of the feedback 
system has been employed and even though the 
method of empirically determining the quan- 
titative characteristics does not give unique 
results. 


BORAX Reactors 


Reference 4 is a report onthe operation of the 
BORAX-IV reactor and covers both the steady- 
state performance and the transient performance 
of the reactor. These experiments gain added 
significance because the reactor has many design 
characteristics in common with the previous 
BORAX reactors and because the effects of the 
specific design changes which were made are 
reasonably well isolated. The fundamental 
change made in the transition from the previous 
reactors to BORAX-IV was the change from a 
highly enriched uranium-aluminum plate type of 
fuel element to an element composed of thoria— 
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Figure 15— Fuel elements of 


highly enriched urania pellets jacketed in alu- 
minum. The additional heat capacity of the 
thoria diluent combined with the relatively low 
thermal conductivity of the thoria-urania sin- 
tered compact resulted in a thermal time con- 
stant which is several times that ofthe previous 
plate type fuel elements. A secondary result of 
the new fuel structure was a steam coefficient 
of reactivity which was rather highly negative 
(comparable to that of BORAX-I), even though 
the BORAX-IV core is relatively large andeven 
though the neutron leakage is relatively low. 
The more negative void coefficient is due pri- 
marily to the high density of U*> in the core, 
which reduces the importance of the neutron 
absorption in H,O, and to the resonance ab- 
sorption in the thorium, which adds a negative 
component to the steam coefficient which was 
not present in the previous BORAX reactors. 
Cross-sectional drawings of the fuel elements 
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the four BORAX reactors. 


used in the series of BORAX reactors (I to IV) 
are shown in relative scale in Fig. 15. 

With regard to reactor stability, the most 
striking result of the experiments was the very 
large amount of excess reactivity which could 
be compensated by steam in stable operation at 
a pressure of 300 psi. For the highest power 
run, this amounted to about 6.9 per cent F off 
and gave a power output of 20.6 Mw(t) from a 
core consisting of 69 fuel elements. The power 
corresponded to a power density of about 69 
kw/liter of coolant. Previous BORAX reactors 
had not been operated with more than 3 per cent 
reactivity compensated by steam, and it ap- 
peared evident in these reactors that further 
additions of reactivity would lead to instability, 
although the largest of the reactors (BORAX-III) 
was never taken to the actual point of instability. 

The power performance of the series of 
BORAX reactors is summarized in Table X-1, 
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Table X-1 CHARACTERISTICS OF FOUR BORAX REACTORS!~* 
A B Cc D E F G H 
Power Power Reactivity D/E 
Operating density in density in compensated F XG, 
pressure, core, Coolant vol. coolant, by steam faved ) Steam coefficient kw/liter 
Reactor psig kw/liter Core vol. kw/liter (AR o6e/R ore) AR os/® ess of reactivity* of steam 
BORAX-I 135 q 11 0.84 13 0.028 465 —0.24F 112 
BORAX-II 300 27 0.702 38.5 0.030 1284 —0.10,¢ —0.19f 128—244 
BORAX-III 300 29 0.711 40.8 0.020 2040 —0.07,¢ —0.10f 143—204 
BORAX-IV 300 54.6 0.056 975 —0.28f 273 
BORAX-IV 300 68.7 0.069 1000 —0.28f 280 








*(AR.¢/R og) /(Steam volume/coolant volume). 
+From reference 5. 
tFrom reference 7. 


which is derived from the data of references 4 
to 7. It is quite evident that very substantial 
increases have been made in the maximum 
power density attainable, as was predicted after 
the earliest boiling-reactor tests. It should be 
noted that the BORAX-I results are not strictly 
comparable with the others, since that reactor 
was operated at 135 psi, whereas the later 
reactors were operated at approximately 300 
psi. If the maximum power density in the coolant 
is divided by the amount of reactivity com- 
pensated by steam at that operating condition, 
the result reflects the variation in steam coef- 
ficient of reactivity from reactor to reactor. 
Thus the quantity increases progressively from 
BORAX-I to BORAX-III, reflecting the decreas- 
ing value of steam coefficient in this series of 
reactors. It was this progressive reduction in 
steam coefficient which accounted for most of 
the improvement in attainable power density in 
the series. In the case of BORAX-IV, however, 
the ratio drops appreciably, illustrating the 
fact that the higher power density of BORAX-IV 
was attained because of the higher permissible 
amount of reactivity compensated by steam and 
was attained despite the higher steam coefficient 
of reactivity. This is equivalent to saying that 
the higher power density resulted at least in 
part from greater reactor stability. 

It appears, however, that the greater stability 
is not responsible for all the improvement in 
performance. If one takes the ratio (kilowatts 
per liter of coolant volume)(steam coefficient of 
reactivity)/(reactivity compensated by steam), 
the result should be the average number of 
kilowatts of heat generated per liter of steam 
contained in the core, and if the rate at which 
steam can flow out of the core were constant, 


one would expect this ratio also to be approxi- 
mately constant. The last column of Table X-1 
shows that this ratio increases gradually over 
the series of BORAX reactors as might be 
expected because of the increasing fuel-plate 
spacing over the series. The implication is that 
an increase in fuel-plate spacing increases the 
velocity of steam flow out of the core. 

It must be stated that too much significance 
cannot be read into comparisons like those 
discussed above in which only the gross varia- 
bles are taken into account. One significant factor 
which is not taken into account and which may 
have a large effect on the stability, as well as 
on the general thermal performance, is that the 
first three reactors were fueled with mixtures 
of fuel elements of rather widely different fuel 
loadings. In some cases this probably led to 
rather high peak-to-average power density 
ratios. Furthermore, the steam coefficients of 
reactivity are calculated values rather than 
measured values and are average values over 
the core rather than true effective values. The 
relation between average and effective values 
would be expected to vary over the series of 
reactors. 

A rather extensive investigation of the sta- 
bility of BORAX-IV at atmospheric pressuresis 
also reported in reference 4. Although this con- 
dition, as an operating condition for power 
reactors, is not an interesting one, it is never- 
theless worth investigating because the stability 
characteristics become poorer as pressure 
is decreased. Several theoretical studies of 
boiling-reactor stability have concluded that the 
stability characteristics of a reactor at atmos- 
pheric pressure are qualitatively different from 
those at high pressure and that nogreat concern 
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need be felt because the low-pressure charac- 
teristics are relatively poor.’ It would, neverthe- 
less, be comforting to understand in detail the 
low-pressure instabilities and to be able to 
predict the operating pressure at which they 
become insignificant and at which other effects 
pecome predominant. 

It appears that, at atmospheric pressure, 
BORAX-IV can tolerate no greater reactivity 
compensation by steam than that tolerated by 
previous BORAX reactors. Indeed, the insta- 
bility which has in the past been referred to as 
“chugging” was observed when the reactivity 
compensated by steam was no more than about 
1.5 or 1.6 per cent, whereas BORAX-I had on 
occasion operated well above 2 per cent at 
atmospheric pressure without encountering in- 
Stability. The power density attainable with a 
given amount of reactivity compensated by 
steam was, however, considerably higher for 
BORAX-IV than for BORAX-I, which was the only 
other reactor for which atmospheric-pressure 
power measurements were made. Reactor IV 
attained a power density of some 25 kw/liter of 
coolant when the reactivity compensated by 
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steam was about 1.2 per cent. With this reac- 
tivity compensation, BORAX-I attained a power 
density of only about 8 kw/liter of coolant. 

The transfer function of BORAX-IV was meas- 
ured at atmospheric pressure over the power 
range from effectively 0 to 3.4 Mw(t). Although 
the measured transfer functions did not have the 
highest precision because of the necessity for 
using relatively crude equipment and because of 
the “noise” inherent in a boiling system, sig- 
nificant data were obtained which are sum- 
marized in Fig. 16. As reactor power was 
increased, a marked variation in the transfer 
function was observed. At powers of approxi- 
mately 1 Mw and higher, a maximum was ob- 
served in the gain curve, corresponding to a 
damped resonance of the reactor power with 
the impressed reactivity oscillation. As power 
increased, the amplitude of the resonance in- 
creased, as did also the resonant frequency. 
The solid line in part A of Fig. 16 is the recip- 
rocal of the gain at the resonant frequency as a 
function of reactor power. The gainis defined as 
the fractional increase in reactor power divided 
by the fractional increase in k.s¢ which produces 
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Figure 16—Behavior of gain and frequency of the transfer function with power for BORAX-IV at 


atmospheric pressure. 
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the power change (gain = (AP/P)/(Akg¢s/R ofp): 
The extrapolation of the line to zero reciprocal 
gain would mark the approximate power at which 
instability would be expected. The dashed line 
in part A of Fig. 16 shows the ratio of the gain 
at zero frequency to the gain at the resonant 
frequency. As reactor power is increased, the 
self-regulating effect of the steam production 
becomes more and more pronounced. That is 
to say, the gain of what may be considered the 
steam feedback loop in the self-regulating sys- 
tem of the reactor increases, and the net gain 
of the system, (AP/P)/( ARe¢s/Repe), decreases 
at the frequencies which are low enough to allow 
the self-regulating effect to operate as astable, 
effective system. Consequently, the zero fre- 
quency gain of the system decreases as power 
is increased, and the ratio of the zero frequency 
gain to the resonant frequency gain decreases 
more rapidly than the reciprocal of the resonant 
frequency gain. Part B of Fig. 16 shows the 
resonant frequency as a function of reactor 
power. The very marked increase in resonant 
frequency with power is something that could 
not be adequately demonstrated by the observa- 
tions which had previously been made on the 
self-excited oscillations of other reactors, anc 
it demonstrates the usefulness of the oscillator 
technique in such studies. It seems probable 
that this change of resonant frequency may be ar 
important manifestation of the basic physical 
processes involved and may make an important 
contribution to the understanding of stability in 
such a reactor. 

Methods of obtaining reactor transfer func- 
tions other than by use of a reactivity oscillator 
were investigated. One of these methods called 
the ringing method is a well-known technique 
in other applications; the method consists simply 
of applying a sudden impulse or perturbation to 
the system (in this case a rapid insertion of 
reactivity) and observing the transient response 
of the system. A transfer function canbe derivec 
by appropriate observations and mathematical 
manipulations of the transient. A second 
method®"? is to derive the transfer function 
from the autocorrelation function for the normal 
fluctuations in the reactor power which are 
caused by the boiling “noise” in the reactor. In 
essence this method is much the same as the 
ringing method, but the perturbations which 
excite the ringing are the boiling fluctuations in 
reactivity, and the autocorrelation analysis is 
necessary in order to separate the response of 


the reactor to a typical impulse of this kind 
from its response to previous and subsequent 
impulses. The methods may be useful for reac- 
tors which are not equipped with suitable reac- 
tivity oscillators. 


Review of ANL Experience 


Reference 7 contains a review of experimental 
information obtained from the boiling reactors 
which have been operated by Argonne National 
Laboratory and also contains a theoretical con- 
sideration of the stability question. The experi- 
mental section of the report includes a useful 
collection of information on the design charac- 
teristics and observed behavior of the EBWR and 
the BORAX reactors, much of which is not 
available elsewhere. This section of the report 
also performs a useful service in redirecting 
attention to some of the dynamic characteristics, 
such as the boiling fluctuations of power andthe 
small-amplitude oscillations which characterize 
the boiling reactor when its power level is well 
below the stability limit. In this connection, 
Table X-2, which summarizes some of the 
phenomena by which one might evaluate the 
degree of stability of a reactor, also gives an 
indication of the characteristic behavior that one 


Table X-2 MEASURES OF EBWR STABILITY’ 





¢ from ringing tests +0.13 
Pp for R = 1.0364 0.34 sec! 
Root-mean-square value of p, i.e., 

vp? (includes noise as well as the 

dominant oscillatory frequency) 0.83 Mw 


Half-width of the resonance in the 
power spectrum from autocorrela- 
tion techniques 

Properties of the transfer function G: 
Value at resonance 
Half-width of resonance 


18 counts/sec 


55.4 db 
0.38 counts/sec 





might observe in the stable operation of a typical 
natural-circulation boiling reactor. The values 
given in the table were those observed for the 
EBWR at a thermal power of 45.6 Mw and an 
operating pressure of 600 psig. 

Some of the symbols of Table X-2 require 
explanation. If a “ringing” test is made on the 
reactor by applying a sudden reactivity increase, 
a power transient of approximately the form 


P(t) = (constant) exp (—¢wt) sin wt 
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may be expected. ¢ is the damping constant for 
the transient oscillation and, as may be seen 
from Table X-2, had a value of 0.13 for the 
conditions in question. When the reactor is 
operating steadily, the power is subject to fluc- 
tuations from the “boiling noise” and, at times, 
from the effects of oscillations. Pp is the 
probability per second that the power will ex- 
ceed a value R times the average power. The 
symbol p is used to designate the power fluc- 
tuation, that is, the instantaneous deviation 
from the average value. Statistical analysis of 
these deviations reported in the reference shows 
that their distribution is Gaussian. The root- 
mean-square value of 0.83 Mw for p is 1.8 per 
cent of the average power. The magnitude ofthe 
reactor transfer function at the resonance fre- 
quency is given as an absolute quantity in 
decibels. The amplitude of the transfer function 
is the ratio of the fractional amplitude of the 
power oscillation to the fractional amplitude of 
the reactivity oscillation, (AP/P)/(ARg¢¢ / Rog), 
and the magnitude in decibels is equal to 20 
logi9(AP/P)/( Aes /Ref¢). This absolute specifi - 
cation of the transfer-function magnitude is not 
universally used. For example, in reference 3, 
the 0-db gain is defined as the amplitude of the 
zero-power transfer function in the region of 
its plateau (at a frequency of approximately 20 
radians/sec in Fig. 13). The absolute amplitude 
of the transfer function in the plateau region is, 
for any reactor, 1/8, where f is the delayed 
neutron fraction. Thus 1/8 is approximately 
140 or about 43 db. The magnitude of the transfer 
function at the resonance point quoted in Table 
X-2, 55.4 db, is therefore about 55.4 — 43 = 12.4 
db above the zero-power plateau value. Thus 
the observed fractional amplitude of the power 
oscillation at the resonance frequency was 
higher by a factor of about 4.2 than the fractional 
amplitude of the zero-power oscillation which 
would have been produced by the same reactivity 
oscillation. 

The theoretical section of reference 7 treats 
the derivation and application of a theoretical 
transfer function which neglects effects caused 
by pressure variations (and therefore is limited 
in applicability to the higher frequencies) 
but which treats the heat-transfer, steam- 
production, and steam-escape processes in 
more detail than was done in reference 3. In 
particular, the boiling and nonboiling regions 
of the fuel channels are treated separately. 
The theoretical transfer functions for the EBWR 


and the BORAX-II, -III, and -IV reactors are 
computed numerically for various operating 
conditions. The EBWR computations agree rea- 
sonably well with the experimental amplitude 
transfer functions. The BORAX-II and -III 
calculations show resonances at about 40 
radians/sec. It is not obvious whether the 
calculation indicates an impending instability 
for these cases, but it can be said that the 
computed frequencies are considerably higher 
than the observed oscillation frequencies just 
before the onset of instability in these two 
reactors. The application of a common stability 
theory to a series of reactors is, however, a 
very worthwhile endeavor, for it willbe difficult 
to have complete confidence in the prediction of 
stability until a consistent theory is available 
which explains the stability limits of all those 
boiling reactors which have been investigated. 
Reference 13 is a theoretical study of the 
dynamics of boiling reactors which treats in 
some detail the dynamics of heat flow, steam 
flow, and the flow of water in the fuel channels 
and in the recirculation path. Reference 13 does 
not consider the effects of pressure changes in 
the steam space above the top reflector of the 
reactor; in fact, the reference is applied pri- 
marily to the understanding of the unpressurized 
BORAX andSPERT experiments. It has generally 
been thought that the hydrodynamic forces would 
be most important in such low-pressure reac- 
tors, and, consequently, such an application of the 
theory is a logical one. The analysis is not 
made by the transfer-function approach, but, 
rather, the analysis sets up the appropriate 
equations for the reactor power and solves 
them by machine computation for specific tran- 
sient operational cases such as the sudden or 
ramp addition of reactivity or the sudden addi- 
tion of cold water. In the case of ramp addition 
of reactivity to the mathematical model of the 
subcooled SPERT reactor, the response of the 
reactor appeared to be smooth, and no insta- 
bility was found up to the limit of the reactivity 
addition which amounted to more than 9 per 
cent Refs. This, of course, is contrary to the 
observed SPERT behavior. It is postulated in 
the reference that the unstable behavior of 
SPERT and BORAX-I at atmospheric pressure 
was due primarily to a hydrodynamic instability 
connected with the natural-circulation boiling, 
which, although it certainly would couple with 
the other dynamic factors in the reactor, could 
produce an oscillation in steam content of the 
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reactor even without such coupling. Such oscilla- 
tion of the steaming rate has been observed in 
laboratory tests on boiling channels but has 
been observed only at steam contents consider- 
ably higher than those thought to exist in the 
BORAX and SPERT reactors at their instability 
points. The reference considers this instability 
from the theoretical standpoint and shows thatit 
becomes a possibility when the steam content 
of the channel becomes high enough that a further 
increase in steam content will lead to a de- 
crease of recirculation velocity, because the 
increment of flow resistance resulting from the 
extra steam content is greater than the corre- 
sponding increment in recirculation head. 

The reference also makes some applications 
of the dynamic theory to a heavy-water boiling 
reactor; the reference also points out some 
possible differences between heavy-water boil- 
ing and light-water boiling reactors. It suggests 
that, whereas the closely packed fuel channels 
of the light-water reactor would be expected to 
act in concert with respect to hydrodynamic in- 
stabilities, then the segregated channels of some 
heavy-water boilers might be expected to act 
individually. 
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FLUID-FUEL REACTORS 





Report of the Fluid Fuel Reactors 
Task Force 


In the fall of 1958, the U. S. Atomic Energy 
Commission organized the Ad Hoc Advisory 
Committee on Reactor Policies and Programs. 
The Committee reviewed the AEC’s civilian 
power program and recommended immediate 
and long-range programs. One of the Commit- 
tee’s recommendations is quoted: “The aqueous 
homogeneous, molten-bismuth, and molten-sal* 
reactors all offer the possibility of reducing 
the cost of the fuel cycle, and the last two offer 
the possibility of high-temperature operation. 
These three concepts for power reactors should 
be critically compared and work concentratec 
on the concept that appears the most promising.” 

In accordance with this recommendation, the 
Evaluation and Planning Branch of the Civilian 
Reactors Office, Division of Reactor Develop- 
ment, organized the Fluid Fuel Reactors Task 
Force in January of 1959. The Task Force hac 
as its assignment and objective the comparison 
of three reactor concepts: the aqueous homo- 
geneous reactor (AHR), the molten-salt reac- 
tor (MSR), and the liquid-metal-fuel reactor 
(LMFR). The comparison was to determine, as 
follows: the present state of development and 
technical feasibility; the technical feasibility of 
breeding; the cost of power for an optimized 
system; and the research and development pro- 
gram necessary to achieve this cost. In order 
to study the concepts, it was necessary to have 
reference designs prepared. This design work 
was done by three project groups, one for each 
concept. The report of the Task Force has re- 
cently become available. ! 


Reference Designs 


1. Molten-salt Reactor. The reference- 
design molten-salt reactor is an INOR-8 vessel 
containing a graphite assembly 12.25 ft in di- 
ameter by 12.25 ft high, through which molten- 
salt fuel flows in vertical channels. The fuel 
salt is a solution composed of 0.3 mole % UF,, 
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13 mole % ThF,, 16 mole % BeF,, and 70.7 mole 
% Li;F. The fuel salt is heated from 1075 to 
1225°F in the core and is circulated from the 
reactor vessel to four primary heat exchangers 
by four fuel pumps. A barren coolant salt is 
used as the intermediate heat-exchange fluid, 
which superheats and reheats steam in a Loeffler 
boiler system. The reactor develops 760 Mw of 
heat. By using 2000-psig, 1000/1000° F steam in 
a reheat cycle, the net electrical output is 318 
Mw(e); the net efficiency is 41.9 per cent. 


2. Liquid-metal-fuel Reactor. The liquid- 
metal-fuel reference-design reactor is fueled 
with a slurry containing approximately 3 wt.% 
ThO,-UO, in bismuth. The reactor core vessel 
is 14 ft in diameter and 14 ft high, containing a 
graphite core and reflector assembly. The fuei 
slurry is heated from 750 to 1050°F througk 
channels in the reactor core and is circulated 
through three primary heat exchangers by three 
primary variable-speed pumps. The reactor de- 
velops 825 Mw of heat. Sodium is used as the 
intermediate heat-transfer coolant. The steam 
from the secondary heat exchangers supplies a 
non-reheat, 2000-psi, 1000°F steam turbine 
which produces 312 Mw of net electricity. The 
net efficiency is 37.8 per cent. 


3. Aqueous Homogeneous Reactor. Three 
concepts of aqueous homogeneous reactors were 
presented to the Task Force: a two-region, 
solution-core, slurry-blanket reactor; a two- 
region, Slurry-core, slurry-blanket reactor; 
and a one-region slurry reactor. Heavy water 
is the fuel carrier and moderator in all cases. 
The design group chose the two-region, solution- 
core, slurry-blanket reactor as the reference 
design. This reactor consists of a 4-ft-diameter 
by 12-ft-long zirconium alloy core tank in an 
8-ft-diameter by 16-ft-long pressure vessel. A 
5 g/liter solution of uranyl sulfate in heavy 
water is circulated through the core of the re- 
actor, and a thoria slurry containing 1000 g of 
thorium per liter is circulated through the blan- 
ket. Heat is removed from the core solution in 
two circulating systems, each of which contains 
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a steam generator and a circulating pump. The 
blanket slurry is recirculated through one simi- 
lar heat-removal circuit. Steam is generated at 
400 psia and 435°F and sent directly to the tur- 
bogenerator. The heat-generation rate for one 
reactor is 380 Mw(t): 320 Mwi(t) in the core and 
60 Mw(t) in the blanket. Three reactors provide 
steam for one turbogenerator in the reference 
station to produce 317 Mw(e) net. The slurry- 
fueled reactors are similar but larger; for each 
type, one 1140-Mw(t) reactor satisfies the re- 
quirements of the reference station. 


Present State of Development 


and Technical Feasibility 


For fluid-fuel reactors, the critical factors 
are fuels and their containers. Table XI-1 lists 
estimated technical difficulties for the three 
systems. The larger ratings represent smaller 
probabilities of technical success. 


Table XI-1 RELATIVE TECHNICAL 
DIFFICULTIES— FUELS AND MATERIALS? 











AHR 
Solution Slurry 
core core 
(reference) (alternate) MSR LMFR 
Materials: 
Internal 5 2 2 1 
External 3 2 1 3 
Total 
materials 8 4 3 4 
Fuel 7 5 2 4 
Total 15 9 5 8 





The internal material for the MSR and LMFR 
is graphite, and the AHR uses Zircaloy. The 
basic container and piping for the AHRis INOR-8 
(79 per cent Ni, 12 per cent Mo, 7 per cent Cr, 
5 per cent Fe), whereas the LMFR employs 
Croloy and the AHR employs stainless steel. 
The relatively poor rating for the AHR on inter- 
nal structure, a value of 5, comes fromthe diffi- 
culty of assuring integrity of the Zircaloy core 
tank for all operating conditions.* The LMFR 
is penalized by mass transfer in the external 
piping loops while stress-corrosion cracking and 
hydrogen embrittlement of stainless steel is a 
problem on the AHR. The low(numerical) rating 





*See the section on the Homogeneous Reactor Test, 
p. 58, for additional details. 


for the MSR on internals and externals is due to 
the good corrosion and mass-transfer proper- 
ties of INOR-8 in the presence of fused salts. 

The MSR has, in addition, a good rating for 
fuel, principally because its fuel is a solution. 
Both the AHR and LMFR have slurry problems, 
the magnitude of which “can hardly be defined 
at this time.’’ The rating of 7, assigned to the 
solution core of the AHR, is caused by the ex- 
treme care necessary to maintain the fuel solu- 
tion within proper chemical control to prevent 
damage to the Zircaloy coretank. The alternate 
AHR, with a slurry core, is ina somewhat better 
situation, since corrosion of Zircaloy in the 
slurry case seems to be substantially less than 
in the case of the solution fuel. The slurry fuel 
for the AHR, however, is considered a more 
difficult problem than that of the LMFR (5 ver- 
sus 4) due to the fact that the aqueous slurries 
have faster settling rates than the liquid-metal 
slurries. On the basis of Table XI-1, the MSR 
has been rated as the “best possibility” for de- 
veloping satisfactory fuel and containing mate- 
rials. 

Maintenance is cited as the most important 
factor influencing the practicability of any of the 
three concepts. The problem basically arises 
from the fact that the fluid-fuel reactor circu- 
lates a highly radioactive stream of fuel around 
the plant. The aqueous homogeneous reactor 
experiments have demonstrated many remote 
maintenance techniques employing water as a 
shield (wet maintenance), but the feasibility of 
dry maintenance has not beenproven. The large 
plants under consideration would consist of 
floors of equipment stacked one above the other; 
successful maintenance would be contingent upon 
remotely operated robots, manipulators, and 
other devices. Another maintenance technique 
would be to spread the plant out on one floor, 
allowing overhead access—a costly procedure. 
{t is stated that: “The liquid-metal, molten-salt, 
and aqueous homogeneous concepts are not 
feasible unless remote welding or substitute 
methods of closure are perfected.’’ Another 
difficult problem, as yet unsolved, is to provide 
remote inspection of the remotely made clo- 
sures. The MSR and LMFR have an advantage 
over the AHR in that the solid fuels, when cool, 
retain most of the fission products. 

A considerable discussion of hazards and 
containment is presented in reference 1. The 
AHR is considered to be potentially the most 
hazardous system because of the large amount 
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of energy stored in the primary system at the 
operating condition. The AHR also produces 
radiolytic gas, an explosion hazard, andis more 
susceptible to nuclear instabilities arising from 
hydrodynamic flow perturbations. The latter is 
caused by the fact that the flow is directed 
through the entire AHR core, whereas the other 
systems employ flow in channels, thereby mini- 
mizing hydrodynamic fluctuations. The ability 
of frozen metallic and salt fuel to retain many 
fission products also is an advantage should an 
accident occur. 

The last area of technical feasibility dis- 
cussed’ is chemical reprocessing of fuel. In 
this area there are some differences of opinion; 
these differences primarily concern the location 
of the processing plant (on-site versus off-site) 
and the amount of process costs. The AHR em- 
ploys the Thorex process, so modified as to 
maintain the desired purity of the heavy water. 
The MSR uses the fluoride volatility process, 
and the LMFR uses a modified Thorex process. 


Technical Feasibility of Breeding 


For comparison three kinds of systems can 
be defined, as follows: 

1. Converters (conversion ratio <0.95) 

2. “Hold own” breeders (conversion ratio ~1) 

3. Doublers (doubling time <20 years) 


The need for these various kinds of breeder re- 
actors depends on how rapidly the production of 
nuclear power expands, on the status of uranium 
reserves available to this country, and on the 
degree of fuel utilization available in other re- 
actors which may be in operation. 

The degree of breeding, as measured by the 
breeding ratio or conversion ratio, depends on 
neutron utilization; of the three concepts, the 
AHR makes the best showing since it can reduce 
parasitic losses to a minimum. The optimum 
breeder reactors were found to be two-region 
machines; a one-region breeder was technically 
feasible but was of too large a size (physical 
and electrical) to be considered very attractive. 
Of the three concepts, only the AHR could be 
considered a doubler, the other two being, at 
best, hold owns in their present forms. The 
possibility of making the MSR and LMFR into 
doublers by reducing the fuel inventory by means 
of internal cooling has been considered. The 
externally cooled reactors are penalized by hav- 
ing large amounts of fuel in the primary circuit 
outside the core. The size of the external cir- 


cuit is affected strongly by considerations of 
remote maintenance and thermal stress. If the 
primary heat exchanger could be incorporated 
in the core structure, the resulting internally 
cooled reactor could have a reduced fuel inven- 
tory and, hence, a shorter doubling time if the 
breeding ratio remained constant. This is not, 
however, an easy development. If it were suc- 
cessful, the resulting reactors might have dou- 
bling times roughly equal to that presently pre- 
dicted for the AHR. 


The above conclusions are predicated on 2 
value of 2.28 for 7 (number of neutrons pro- 
duced per neutron absorbed) of U’**. The feasi- 
bility of thermal breeding is affected strongly 
by the value of 7, not only at thermal but at 
epithermal neutron energies, since the MSR and 
LMFR have substantial amounts of nonthermal 
fission. British measurements in DIMPLE and 
GLEEP have indicated considerably lower values 
of 7 for U***, (For a more detailed discussion 
see Power Reactor Technology, 2(2): 16 (March 
1959).) However, recent U. S. results derived 
from critical experiments’ support the thermal 
value 2.28 or a slightly higher one. Reference 1 
concludes that all the fluid-fuel systems could 
probably at least be hold own breeders, even 
with the uncertainty in 7, since, at an 7 value of 
2.28, conversion ratios were 1.05 for the LMFR 
and MSR and 1.07 to 1.10 for the AHR. 


Cost of Power 


The cost studies were done to compare the 
relative cost potential of the three reactor sys- 
tems, and it is stated that the cost data should 
not be compared with costs of conventional or 
other nuclear plants, since the emphasis in the 
studies was placed on relative, rather than abso- 
lute, numbers. Table XI-2 lists the power costs 
presented by the Task Force. 


As indicated previously, the design groups 
disagree with the Task Force on chemical proc- 
essing costs for the MSRandLMFR. The LMFR 
project estimates that processing can be done 
for “much less than the Task Force estimates,”’ 
and the MSR project estimates are about 1 mill 
lower on processing costs. Although, if correct, 
these would place the nonaqueous reactors in an 
even more advantageous position costwise, the 
Task Force states that the differences in power- 
cost estimates among the three concepts are 
so small in relation to the probable accuracy of 
the estimates that they are not considered ade- 
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Table XI-2 SUMMARY OF POWER COSTS 
(FOR REFERENCE REACTORS: 333 Mw(e), GROSS) 











Mills/kw-hr 
LMFR MSR AHR 

Power-plant investment §.72 6.24 6.70 
Chemical- and waste-disposal- 

plant investment - 0.71 0.71 0.89 
Fuel inventory use and burn-up 1.36 1.37 0.90 
Chemical-plant operation and 

maintenance 0.76 0.76 0.96 
Power-plant operation and 

maintenance 1.46 1.57 1.66 

Total power costs, gross 10.0 10.7 BL.I 

Total power costs, net 10.7 y 5 Re 11.5 





termining factor. It should be pointed out that 
the AHR reference design is for a breeder, 
whereas the other two reactors are converters. 


Research and Development Programs 


The time and cost estimates for research and 
development programs are given in Table XI-3. 


Table XI-3 TIME 


57 


program, to be known as the Thermal Breeder Re- 
actor Program, has been assigned to the Commis- 
sion’s Oak Ridge Operations Office. 

The new program has as its objective the devel- 
opment of a thermal breeder reactor capable of 
converting thorium to fissionable fuel material at a 
doubling time of not more than 25 years. 

The new program also involves a reorientation of 
work on the Fluid Fuel Reactor Program previously 
conducted for the Commission by the Brookhaven 
National Laboratory, Upton, N. Y., and The Bab- 
cock & Wilcox Company, Lynchburg, Va., on the 
Liquid Metal Fuel Reactor (LMFR) concept and by 
the Oak Ridge National Laboratory, Oak Ridge, 
Tenn., on the Aqueous Homogeneous Reactor (AHR) 
and Molten Salt Reactor (MSR) concepts. Much of 
the work on these projects has been directed toward 
the development of reactors of these concepts for 
the production of economic electric power. During 
Fiscal Year 1960, research and development on the 
AHR, MSR, and LMFR concepts, except as applica- 
ble to the new program, will be discontinued and a 
transition of financial support made to the new pro- 
gram. 

The Commission’s decision to reorient the Fluid 
Fuel Program and toward implementation of the 


AND COST ESTIMATES 














AHR MSR LMFR 
Cost, Cost, Cost, 
108 108 108 
Start End dollars Start End dollars Start End dollars 
R&D for reactor experiment July 18 July 19 April 14 
1963 1963 1962 
Design, construction, and July July 15 July December 18 September 14 
operation of reactor 1961 1966 1959 1965. 1964 
experiment 
R&D for prototype 1962 1967 35 1963 1968 35 1964 17 
Design, construction, and 1964 1969 59 1965 50 1963 1971 
operation of prototype 45 
R&D for full-scale plant 1968 1971 5 1968 1971 Not 1965 1969 10+ 
given 
Total 132 122(+?) 100+ 
Conclusions Thermal Breeder Reactor Program was based on 


With the completion of the Task Force study, 
evaluation of the fluid-fuel concept, as compared 
with other power reactor concepts, should be 
possible. The following quotations, taken from 
an AEC press release, outline a newprogram to 
develop thermal breeder reactors:° 


‘‘The Atomic Energy Commission has initiated a 
long-range program to develop effective thermal 
breeder reactors which would make full use of the 
latent energy in thorium. Responsibility for the 


the results of a study by a special task force of 
scientists and engineers convened by the Commis- 
sion in January 1959 to make a comparison of the 
reactor concepts included in the Fluid Fuel Pro- 
gram. As a result of task force estimates of the 
power costs of large-scale fluid-fuel reactors and 
on the cost and scope of research and development 
compared with that of other reactor concepts in ad- 
vanced stages of development, the Commission con- 
cluded that there is little incentive for further de- 
velopment of fluid-fuel reactors as an approach to 
the attainment of low-cost nuclear power in the near 
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future. Therefore the Commission decided to re- 
orient its fluid-fuel effort toward the use of the la- 
tent energy in thorium and natural uranium through 
a long-range development program. Because all 
three of the fluid-fuel reactor systems were found 
by the task force to be potentially capable of de- 
velopment as breeder reactors, the new Thermal 
Breeder Program may include one or more of the 
concepts included in the Fluid Fuel Program. 

In the initial phases, the new program will be a 
basic research and development effort and will in- 
clude an evaluation of reactor technology to deter- 
mine the most promising type of thermal breeder 
reactor for continued long-range development. Any 
work continued on present fluid-fuel reactor con- 
cepts will be conducted in the framework of the 
broader Thermal Breeder Program. 


The Homogeneous Reactor Test 


The AHR concept was reviewed in a previous 
edition of this publication.’ Since that review, 
additional technical information has been accu- 
mulated pertinent to the reactor type, the data 
coming from the Homogeneous Reactor Test 
(HRT) operations and supporting development 
work.*-* The previously discussed Task Force 
report cited the materials and fuel stability 
oroblems of the AHR as being major problems, 
and therefore this discussion will be oriented 
toward those areas. 

The first power run of the HRT, run 13, was 
started on Feb. 7, 1958. This run was initiated 
without any copper recombination catalyst inthe 
fuel so that recombination rate constants could 
be evaluated as copper was added. High 
stainless-steel corrosion rates were experi- 
enced (up to 166 mils/year), but these were re- 
duced to about 1 mil/year in subsequent runs by 
the addition of copper sulfate and an increase in 
the oxygen feed to the stream. Run 13 was ter- 
minated on February 13 when the nickel con- 
centration became too high for reliable chemical 
analysis. Run 14 was initiated with fresh fuel 
on March 24; during the intervening month, 
about 1000 ft of electrical wiring was replaced 
due to faulty insulation. The wet maintenance 
philosophy made this possible without a single 
case of personnel overexposure. On April 4, 
the power was increased to the nominal design 
value of 5 Mw(t). After 0.5 hr at this power a 
leak developed between the core and blanket 
regions of the reactor through the Zircaloy core 
tank, and the reactor was shut down. The leak 
was subsequently determined to be a nearly 


round hole between 1 and 2 in. in diameter at 
the elevation of the weldment of the 30- and 90- 
deg flow transition cones at the bottom of the 
core tank. Mounted in these transition cones, 
perpendicular to the flow direction, were a total 
of nine Zircaloy flow-diffuser screens. Probing 
revealed that the first, second, and fifth screens 
were not attached to the core. Runs 15 and 16 
were undertaken to establish the feasibility of 
operating the reactor in the damaged condition 
with fuel in both the core and blanket regions. 
These runs proved that such operation was 
feasible. 

During the several runs, power surges had 
been experienced; the most drastic run was 
from a power set point value of 3.8 to 5.6 Mw. 
In addition, chemical analysis of the fuel stream 
showed that the concentration of uranium in cir- 
culation was erratic and did not correspond to 
the amount that should have been present by 
material balance. The amount of “lost” uranium 
appeared to increase with increasing power. 
Run 17 was undertaken to examine these prob- 
lems further and to expose corrosion samples 
that were placed in the core and blanket regions 
prior to the run. 

At the end of the run it was ascertained that 
22 of the original 48 corrosion specimens had 
become separated from the titanium sample 
holder and had fallen onto the top diffuser 
screen. Severe corrosion was observed on the 
top screen and adjacent core region, and a hole 
cutter was used to cut 2-in. samples out of the 
top six screens. It was concluded’ that: 


1. Localized sections of the titanium specimen 
holder had been heated to a temperature above. the 
melting point during reactor operation, causing a 
reaction between the holder and the steam—heavy 
water environment. 

2. When the heat source (probably a deposit con- 
taining a high concentration of uranium) was re- 
moved as a result of the reaction, the reaction 
quenched out, indicating that the reaction was not 
self-sustaining. 

3. The damage observed on the Zircaloy-2 
mechanical-property specimens, originally sup- 
ported on the titanium specimen holder, was caused 
by the **burning’’ of the specimen holder and not by 
a reaction of the specimens with their environment. 

4, The damage observed on the top screen of the 
core tank, due to a high-temperature excursion, 
was caused by anevent associated with the titanium 
specimen-holder spike. 

5. The exfoliation observed on the screen speci- 
mens was caused primarily by corrosion or hot 
spots occurring in the cracks and laminations that 
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existed in the screens as installed in the vessel. 
The cracks and laminations were formed primarily 
during the punching of the holes during manufacture. 


With respect to the fuel stability problems, the 
following data® are applicable. 


It is concluded that uranium separates reversibly 
from the fuel solution in a manner which is related 
to the composition of the.fuel and to the reactor 
power. From the criticality data, it appears that 
much of the separated uranium collects on the metal 
surfaces of the core tank, presumably because they 
are at high temperature. After the fuel solution has 
lost a kilogram or more of uranium, power excur- 
sions occur that can be attributed to shifting of the 
material in the core or to redissolution with result- 
ing changes in reactivity. If enough uranium col- 
lects in one place on a metal surface in the core, 
the nuclear heating causes the temperature to rise 
uncontrollably, and the metal melts. Examination of 
the interior of the core with periscopes and metal- 
lographic examination of the corrosion specimens 
and pieces removed from the perforated plates in 
the core inlet indicate that some of the metal has 
been heated to temperatures above 1000°C. There 
is reason to believe that the hole in the core tank 
and the severe corrosion of the perforated plates in 
the core inlet resulted from overheating. 


Run 18 was designed to study fuel instability 
problems and was conducted at a maximum 
power of 1.6 Mw. No power excursions were de- 
tected above background, about 400 gofuranium 
was “lost,” and changes in fuel composition 
were small. 

Run 19 was started to investigate core and 
blanket circulation rate, but low flow in both 
systems caused shutdown. It was found that the 
core circulation pump had been electrically 
connected in reverse. Examination of the blanket 
circulation pump revealed that about three- 
fourths of the titanium impeller had been con- 
sumed, apparently by combustion. 


Although the exact circumstances which resulted 
in the damaged titanium probably cannot be deter- 
mined, it is conjectured that the pump became gas 
bound in a high oxygen atmosphere, and that the 
protective oxide layer onthe titanium was scratched, 
exposing metal to the oxygen. Newly exposed tita- 
nium metal can be expected to burn in an oxygen 
atmosphere.° 


Several designs for areplacement for the core 
and pressure vessel have been prepared. They 
provide for the sweeping of the inside of the core 
tank with inlet fuel in one design andalso for the 
using of the blanket as a core tank coolant in 


another design. Using the blanket as acore tank 
coolant would be accomplished by fabricating the 
core tank with two concentric spherical shells 
and then by pumping blanket material through the 
annular space before discharging the material to 
the main blanket region. In the interim, consid- 
eration is being given to the repair ofthe exist- 
ing reactor by plugging the hole in the core tank 
as follows:° 


1. Thescreen system will be surveyed with peri- 
scopes to determine which screens are already 
loose and to note any special problems. 

2. The screens will be cut loose from the wall. 

3. By holding the loose screens out of the way, 
photographs and acast impression of the hole in the 
core tank will be made to enable a well-fitting patch 
to be machined. 

4. The loose screens will be sectioned or dis- 
solved and removed. 

5. The hole will be plugged. 


Needless to say, this procedure would require a 
considerable number of special tools, such as 
remote cutters, millers, nibblers, and loose- 
screen manipulators. A special core-hole im- 
pression device has been designed to press a 
thermoplastic material into the hole, to heat the 
plastic to soften it, and then to cool the mold 
prior to removal from the hole. Another tool 
has been constructed to toggle-bolt a patch over 
the hole using a gold gasket for sealing. The tool 
has been successfully used in the HRT mock-up. 
An ultrasonic thickness gauge to measure core- 
wall thickness has been built and tested suc- 
cessfully. Another interesting device is the re- 
mote weld-impression device, designed to take 
an impression of the welds in the top of the 
pressure vessel. This is accomplished by em- 
ploying a block of thermoplastic which is heated 
by an embedded resistance heater. 

If a successful repair is accomplished, it is 
obviously of more than passing interest to as- 
sure that another hole is not burned in the core 
tank. Hydraulic tests in the HRT mock-up have 
resulted in the following conclusions:°® 


As a result of the above investigations, it was 
concluded that satisfactory hydrodynamic operation 
of the HRT core could be achieved only by removing 
the top five or more screens and reversing the di- 
rection of flow. The importance of preventing solids 
deposition in a region of high power density may be 
appreciated by considering that the power density 
in corrosion products containing 4 per cent uranium 
is approximately 20 times that of the fuel solution. 
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NEW REACTOR CONCEPTS: THE FLUIDIZED-BED REACTOR STUDY 





The fluidized-bed reactor (FBR) has been under 
study at Martin Nuclear Division for some 
time, and reference 1 is a report of recent 
progress. The reactor primarily consists of a 
cylindrical tank with retaining screens placed 
across each end: When operating the reactor, 
fuel pellets are placed in the tank and water 
is pumped into the tank from the bottom. At 
a certain velocity the drag of the water (or 
other coolant) flowing through the bed com- 
pensates gravitational attraction, and the pel- 
lets become suspended; also, they are con- 
tinuously agitated by the turbulent flow of 
coolant. Reference 1 contains fluid-flow, heat- 
transfer, and nuclear analyses directed at 
determining the feasibility of such a concept. 
The preliminary studies were so sufficiently 
promising that a reference reactor was de- 
signed. Pertinent design data are presented in 
Table XII-1. 

The pressure vessel is constructed of SA-353 
clad with 304 stainless steel. A bottom grid 
retains the pellets upon shutdown, and a top 
retaining grid prevents loss of fuel in the event 
of pump runaway. A particle separator is in- 
corporated into the structure above the core; 
this is a cyclone separator which helps to 
remove fuel particles from the coolant before 
it leaves the vessel. Provisions are made for 
continuously slurring pellets from the core for 
chemical processing and for returning fresh 
fuel pellets to the core. 

The control-system analysis for this reactor 
is particularly interesting since the nuclear 
behavior is closely related to the flow rate of 
coolant through the parameter of bed porosity. 
A considerable number of analog-computer 
runs were made to ascertain whether the re- 
actor could be shimmed and/or scrammed by 
controlling the primary-coolant flow rate. The 
scram conditions analyzed were a cold-water 
accident, a high-power scram trip event, and 
a ramp insertion of positive reactivity. The 
response of the system, using flow rate con- 
trol only, to increase power from 35 to 90 per 
cent, was determined. Within the assumptions 


Table XII-1 FBR DESIGN DATA! 





Reactor power 


Gross electrical output, Mw 

Steam pressure at full load at steam 
generator, psia 

Normal reactor—coolant—system operating 
pressure, psia 


1366 x 10° 
Btu /hr; 
400 Mw(t) 

117 

615 


2000 


Reactor—coolant-system design pressure, 2500 
psia 
Total number of loops 3 
Steam generators per loop 1 
Pumps per loop 1 
No. of loops operating for design power 3 
of 117 Mw 
Flow per loop at design power, Ib/hr 
Estimated reactor pressure drop, psi 75 
Estimated steam-generator pressure drop, 45 
psi 
Estimated primary-pipe pressure drop, psi 20 


Estimated throttling, psi 40 
Total estimated pressure drop, psi 180 
Steam-generator inlet temperature, °F 365.5 
Steam-generator outlet temperature, °F 491.5 
Velocity in main loop piping at 528.0°F, 30.2 
ft/sec 
Reactor coolant flow (total), Ib/hr 15.21 x 10° 
Reactor inlet temperature, °F 491.5 
Reactor outlet temperature, °F 365.5 
Steam-generator size (total), sq ft 55,500 
Primary-piping diameter, in. 16 
Pressure-vessel diameter, ft 105 
Pressure-vessel length, ft 13 
Weight of fuel, Ib 63,736 
Core diameter, ft 8.5 
Core height, ft 6.0 
Fuel—pellet material UO, 
Fuel-pellet diameter, mils 250 
Enrichment, % 15 
Volume fraction of coolant 0.70 





made in solving the system-response equations, 
the technique of controlling the reactor by 
controlling the primary-coolant flow appeared 
satisfactory. 

The reactor contains about 64,000 lb of UO,, 
with a resulting specific power of about 12.5 
Mw(t) per ton of oxide. This is typical of other 
water-cooled reactors of lowenrichment (Dres- 
den: 10.5, Yankee: 21.6). The specific power 
attests to the good heat-transfer properties of 
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fluidized beds. Although the superficial coolant 
velocity* was about 1.40 ft/sec in the bed, the 
large heat-transfer area and large amount of 
turbulence promoted high heat-transfer rates. 
One of the main advantages sought in this 
reactor type is reduced fuel cost primarily 
through reductions in both fabrication and proc- 
essing costs; however, improvements of neu- 
tron economy, which would reduce the net 
burn-up cost, are also visualized. The amount 
of cost improvement which these anticipated 
reductions represent depends strongly on the 
attainable fuel life in the more conventional 
reactors. The Yankee reactor fuel costs, for 
example, are composed of roughly 40 per cent 
element fabrication costs and 60 per cent proc- 
essing, burn-up, and transportation costs. Of 
these costs, both the fabrication and processing 
costs, when referred to unit energy production, 
vary approximately inversely with fuel life- 
time. These variations are illustrated in Table 
XII-2, which shows the estimated Yankee fuel- 
cost breakdown as a function of fuel lifetime. 


pellets caused by the relatively violent con- 
ditions in the bed. Erosion and corrosion of 
UO, pellets by water alone was found to be in- 
significant at the operating conditions of the 
reactor. The problem, of course, can be partially 
solved by minimizing the effect of fines in the 
system as well as by securing a more abrasion- 
resistant pellet. The former could be accom- 
plished by processing a side stream for re- 
moval of fines, designing the primary loops to 
ensure no “dead spots” for trapping of fines, 
providing efficient fines removal equipment in 
the reactor cyclone system, achieving satis- 
factory decontamination techniques, etc. The 
current design of the FBR called for a maxi- 
mum of 500 ppm fines in the main coolant 
loops. Preliminary abrasion tests on UO, pel- 
lets in small fluidized beds at room tempera- 
ture indicated that the pellets lost from about 
0.5 to 1.5 per cent of their weight per day. A 
1 per cent loss yields a fines concentration of 
4000 ppm in the reference design, and hence 
abrasion resistance must be improved if the 


Table XII-2_ YANKEE FUEL-CYCLE COSTS? 





Cost, mills/kw-hr 








Chemical 
Irradiation Initial Fuel- processing 
level, ue element and uranium p35 Pu 

Mwd/ton assay, % fabrication Transport conversion burn-up credit Total 
2,000 2.33 7.44 0.74 2.01 3.28 —-1.32 12.15 
4,000 2.48 3.72 0.37 1.00 2.64 —1.14 6.59 
6,000 2.60 2.48 0.25 0.67 2.50 —1.08 4.82 
8,000 2.77 1.86 0.19 0.50 2.36 —-1.01 3.90 
10,000 2.90 1.49 0.15 0.40 2.29 —0.94 3.39 
12,000 3.05 1.24 0.12 0.33 2.23 —0.85 3.07 
14,000 3.20 1.06 0.11 0.29 2.17 —0.77 2.86 








A key point in the fluidized-bed reactor con- 
cept is the fuel element. The economic success 
of the reactor could well depend on obtaining 
a fuel pellet that did not seriously contaminate 
the coolant system through the release of fission 
products by one means or another. The im- 
portance of obtaining satisfactory fuel pellets 
was recognized in the feasibility study, and 
considerable space was devoted to a discussion 
of fuel pellets. A major problem connected with 
the unclad UO, pellet fuel is abrasion of the 





*The superficial velocity is the velocity the fluid 
would have, at a given mass flow rate, if there were 
no particles in the container. 


design condition is to be met. It was recognized 
that testing at elevated temperatures and in 
radiation fields is a necessary part of a fuel- 
pellet development program. 

Several experiments to clad the UO, pellets 
with nickel were described. The process is 
called “electroless nickel plating” and is ac- 
complished as follows: 

1. The UO, pellets are cleaned in hot alkaline 
solution and washed in distilled water. 

2. Dilute palladium chloride is used to acti- 
vate the UO, surface. 

3. The pellets are placed in nickel sulfate — 
sodium acetate—sodium hypophosphate for 1.5 
hr. 
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The resulting nickel cladding, of about 1 mil 
in thickness, is mechanically bonded to the 
UO, surface. Although the cycling and abrasion 
tests were of fairly short duration, several 
thermal-cycling tests and abrasion tests were 
reported, and the results were encouraging. 

The UO, pellets were prepared by two meth- 
ods. The first method:was one of cold com- 
pacting followed by sintering, and the second 
method was a capillary-drop method. This 
latter method consists of dropping a water 
suspension of UO, into a water repellent bed 
and drying the drop while it retains a spherical 
shape; these events are followed by sintering. 
Both methods produced UO, pellets having den- 
sities of about 90 per cent of the theoretical 
density. Additives were tried in both methods 
in an attempt to improve the physical prop- 
erties. The additives studied were TiO,, CaO, 
CaF,, and CeO,; these additives increased 
sintered densities to about 95 per cent. It 
was also hoped that the additives would permit 
increased bonding at the grain boundaries lead- 
ing to increased abrasion resistance. ThO,-UO, 
pellets were also made by the cold-compacting 
process and had densities of 98 per cent of 
the theoretical density following sintering. 

As pointed out in reference 1, either jacketed 
or unjacketed fuel pellets might be considered 
for a fluidized-bed reactor, and much of the 
investigation summarized in the reference cov- 
ers problems which would be common to both 
systems. However, the two systems would differ 
greatly in many of their characteristics and 
problems; the systems should probably be con- 
sidered as two distinct reactor types. 

If jacketed fuel elements were used, the pos- 
sible advantages of the fluidized-bed system 
would include the elimination of control rods, 
the compensation of fuel burn-up by the loading 
of new fuel in small increments, some re- 
duction in parasitic neutron absorption, and 
perhaps lower unit cost for fuel fabrication. 
An inherent problem, in addition to those prob- 
lems (which were investigated in the reference) 
stemming directly from the characteristics of 
a fluidized bed, would seem to be that there 
is no obvious method of selecting old or dam- 
aged fuel pellets for removal as new fuel is 
added. Apparently, in such a reactor one would 
approach an equilibrium situation in which new 
fuel pellets are added at appropriate intervals 
and fuel pellets are removed in whatever order 
may be dictated by chance. It would appear that, 


in such a system, operation with some fraction 
of defective fuel elements would be necessary; 
that the average fuel lifetime, whether deter- 
mined by fuel-pellet failure or by burn-up of 
fissionable material, would be considerably 
shorter than the maximum lifetime of individual 
pellets. 

If, however, the use of unjacketed ceramic 
fuel in a water-cooled reactor is deemed at- 
tractive, then the fluidized-bed approach is 
one which naturally suggests itself, since the 
more conventional arrangement of fuel in a 
fixed solid structure is probably not easy to 
attain without the use of a metallic supporting 
structure of some kind. The problems of fission- 
product contamination which accompany unjack- 
eted pellets must then be faced. The problem 
of selecting fuel for removal may be automati- 
cally solved, if attrition of the fuel by abrasion 
is severe; presumably, the fuel, on the average, 
will wear away at a more or less constant rate. 
Indeed, if the abrasion rate amounted to ap- 
proximately 1 per cent per day, as in the case 
of the pellets investigated in the reference, 
the average fuel life in the reference reactor 
would be only about 1200 Mwd/ton, even if fuel 
pellets could be left in until they were com- 
pletely worn away. It is doubtful that such a 
short fuel life could be tolerated, even if fuel 
fabrication costs were quite low. Consequently, 
the improvement of abrasion resistance ap- 
pears to be necessary from the standpoint of 
fuel cost as well as from the standpoint of 
system contamination. 

The discussion of pellet fuel elements for the 
Pebble Bed Reactor’ (reviewed in Power Re- 
actor Technology, 2(3): (June 1959)) may also 
be pertinent to the FBR, since there appears to 
be no fundamental reason why UO, is the only 
material that can be used for preparing fuel 
pellets. 
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